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XCVII. An Unbranched Laminar Model of the Intermediate State of 
Superconductors. 


By C. G. Kurpsr, 
Royal Society Mond Laboratory, Cambridge*. 


[Received March 20, 1951.] 


ABSTRACT. 


The experiments of Meshkovsky and Shalnikov on the domain structure 
of a sphere in the intermediate state give no support to the “ branching ” 
model described by Landau. The present theory assumes an unbranched 
laminar structure, and does not involve any precise assumptions about the 
shape of superconducting regions near the surface of the specimen. 
Assuming only that demagnetizing coefficients for the superconducting 
domains can be defined, the Landau thermodynamical potential 
L=F—(HB+-hB—hH)/87 is found. By minimizing this function, the 
spacing and size of the superconducting domains can be calculated. The 
magnetization curve for ellipsoidal specimens is calculated, and is in 
moderate agreement with experiment. 


§1. IvrRopvucTIoN. 


WHEN a superconducting body with non-zero demagnetizing coefficient 1 
is placed in a magnetic field h>(1—n)H¢ (where Hg is the critical field for 
a long rod) the body must split up in some way into normal and super- 
conducting regions. For a body of macroscopic dimensions these regions 
must be small in comparison with the body, as no large-scale equilibrium 
configuration exists. However, the structure is not indefinitely fine, as 
there may be a free energy « per unit area associated with an interface 
between normal (N) and superconducting (S) regions (London 1937, 
Peierls 1936). It is convenient to work in terms of a quantity of the 
dimensions of a length 4=87«/H72,, rather than in terms of « itself. 

In 1937, Landau put forward the first detailed theory of the structure : a 
laminar model with the S-layers rounded at the free surface of the specimen 
(fig. 1a). On the interface between any N- and S-laminz he assumed 
H=H,. But later he showed that difficulties arise from the curvature of 
an interface and, to avoid these, proposed a new model (Landau 1938, 
1943), in which the N-regions progressively subdivide as they approach 
the surface, so that the domain pattern does not appear at the surface 
(fig. 16). Recent experiments (in which the actual magnetic field in a 
diametral plane and at the surface of a sphere was plotted by traversing it 
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with a small bismuth probe) show, however, that an elaborate domain 
pattern does appear at the surface (Shalnikov 1945, Meshkovsky and 
Shalnikov 1947a, b, Meshkovsky 1949). | Meshkovsky does find occasional 
branching, but it is so rare that it is probably of little quantitative 
importance. 

Andrew (1948 b) gave a theoretical calculation of the magnetization 
curves of superconducting ellipsoids in the intermediate state, both on 
the basis of Landau’s branched Jaminar model and on a similar branched 
thread structure. Comparison with the experiments of Désirant and 
Shoenberg (1948) on thin cylinders and of Andrew and Lock (1950) on 
thin plates, in transverse fields, gave moderate agreement but with 
some noteworthy divergencies. Owing to difficulties in interpretation, 
the resistance measurements of Andrew (1948a) give rather less 


Fig. 1. 


(a) (b) 
Landau’s unbranched Landau’s branched model. 
laminar model. 


information than the magnetization curves, but appear to accord reasonably 
well with the magnetization measurements of Désirant and Shoenberg for 
the same specimens. In view of the above experimental evidence that 
branching does not occur, similar calculations using an unbranched 
model, more like Landau’s original (1937) structure, have been carried out. 

The objections raised by Landau (1943) to his unbranched model (1937) 
appear less serious than when first raised. For if the interface moves into 
the region (A, fig. 1 (a)) where the field is less than H,, while it is true that 
the free energy in the specimen is lowered, the field energy is raised, and so 
we should expect to find an equilibrium configuration. 

In Landau’s (1937) theory, despite its mathematical complexity, the 
simplifying assumption was made that H=H, on an N-S interface. 
However, it can readily be shown (as in Kelvin’s treatment of vapour 
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pressure over curved surfaces) that H¢/H,=1—A/2r, where Hi, is the 
modified critical field and 1/r is the Gaussian curvature. A correction on 
these lines ought also to be made for the sharp corners involved in Landau’s 
branching model, and for the penetration of the field into the thinner 
regions of the superconducting domains. 

'The present theory, on the other hand, does not involve a detailed 
assumption about the shape of the interface near the free surface. It deals 
instead with the demagnetizing coefficients and the surface/volume 
ratios of S-domains. 

In §2 we begin by writing down the Landau thermodynamic potential 
(Landau 1943), which is shown in the appendix to be the appropriate one 
if h, the applied external field, is specified, 


L=Jdt{F—(#HB+NB—-hH)|8r, 


for an unbranched laminar structure. We then minimize L with respect 
both to the number of S-domains and their size. The condition for the 
minimum will be used in §3 to compute the magnetic behaviour of 
ellipsoidal specimens. The Landau thermodynamic potential for the 
present model has a lower value than for Landau’s “ branched ”’ theory 
over a part of the range of applied magnetic field. Over the remainder of 
the range the unbranched state may perhaps be metastable, but this 
question is not investigated here. The above procedure restricts us to an 
equilibrium theory (as for Landau’s 1943 theory). However, the experi- 
ments of Andrew and Lock (1950) show that irreversible effects may be of 
considerable importance. After investigating the equilibrium properties 
we shall discuss in what circumstances hysteresis is to be expected. 

We need not make any explicit allowance for penetration of magnetic 
flux into the S-regions, for it has been shown that, to the first order in 
A/r, the effect of penetration is merely to introduce an additional surface 
energy AHZ2/87. (A is here the penetration depth (Ginsburg 1945, 
Pippard 1951.) If y is defined as the “ true ” total surface energy per 
unit area when allowance is made for penetration, then 4=(87y/H?)—). 
Our A corresponds to Ginsburg’s «’,,,. 87/HZ, or Désirant and Shoenberg’s 4’. 

As in the older theories, we will only consider models in which the N-S 
interfaces are parallel to the applied field except near the surface of the 
specimen. The justification for this is that any other configuration 
would increase the field energy without any compensating reduction in 
the interfacial surface energy. In the plane perpendicular to the applied 
field the structure is found experimentally to be complex—both thread- 
like and warped laminar domains seem to be present (Meshkovsky 1949). 
This suggests that the Landau function is rather insensitive to details of 
the structure, a conclusion borne out by Andrew’s (1948 b) comparison 
of branched thread and laminar models. Any model amenable to calcula- 
tion must be an idealization, and we shall choose a laminar structure for 
simplicity. 

auz2 
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§ 2. THE UNBRANCHED LAMINAR MODEL. 
We shall use the following notation : 


h=the applied magnetic field. 
H =the local value of the magnetic intensity at any point in the body 
B=the local value of the magnetic induction. 
H=the mean magnetic intensity, taken over a region large enough 
to include several domains. 
B=the mean magnetic induction, taken over a similar region. 


We use the macroscopic description of a superconductor in which the 
magnetic fields of the persistent currents are allowed for by assuming a 
difference between Z and # in the S-domains, rather than by the London 
(1937) method of specifying the currents themselves and taking B= H 
everywhere. The two pictures are equivalent provided the S-regions are 
simply connected. 

Temperley (see appendix) has investigated the use of the Landau 
thermodynamic potential and has shown that it is expressible in the form 


Heder | N-regions He dr | interfaces AH dS , (1) 


Provided we are dealing with a body fairly large compared with the 
domain size, it is unnecessary to make detailed assumptions about the form 
of the interphase surfaces near the edge. The only actual assumption 
is that the domains are sufficiently nearly ellipsoidal to be described in 
terms of a demagnetizing coefficient. For this purpose we may regard the 
regions outside the body where the field is very small (A, fig. 2) as part of the 
appropriate S-domains, but if the width of the domain is small compared 
with the thickness of the specimen, this correction will be of little 
importance. 

The method to be used will apply only to ellipsoidal specimens, and to 
them only on the assumption that the density of S-domains is homogeneous 
throughout the body. If this assumption is true, then H is constant 
throughout the body and the body can be described in terms of 
a demagnetizing coefficient m, such that 


87L—const.+h | 


S-regions 


/ 


h=nB+(1—n)H . SM tk IA fool b): 


If A is the cross-section of an S-domain and p is the number of S-domains 
per unit area of cross-section of the specimen, 


B=H(1—pA). Rae ue) fC). 


Ifvis the demagnetizing coefficient of an S-domain, and # is the magnetic 
field within it, 


Hes (1—v) 3 SO Re eile ee mete 
B is zero since we are dealing with an S-domain. Hence 


HHI v)(1—npAyy. a 4 EY 
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Then the Landau function per unit volume, apart from constant terms, is 
given by 


849 L=L=h*pAl|(l—v)(1—npA)]|—H?2,pA(1—As/A) . . . (6) 
or, defining 7»=h/H , 
L = HopA {y?/[(1—v)(1—npA)]—( (L=As/A) 5 5 (6) 


S$ is the perimeter of an S-region in the cross-section perpendicular to 
the applied field h- 


sek 
Fig. 2. 


Field distribution for an unbranched model, showing the external 
regions (A) where the field is small. 


Ginsburg’s (1945) surface energies at the free surface of the body have 
been neglected. According to recent work by Pippard (1951) the difference 
between the N and S surface energies is in any case zero. 

Minimizing Y with respect to p, we find 


(1—npA)?=r?/[(1— 4s/A)(1—v)] 
or npA=1—y/4/[(1— 4s/A)(1—v)]. cane () 
Substituting (7) in (6) and simplifying gives 
—(1/n)[h[x/(1—v)-Hop/(I— As/A)}2. ss (8) 


For further progress it is necessary to specialize the model. We shall 
discuss the case of lamin, since the expression for the demagnetizing 
coefficient for a thread is mathematically more complicated, involving a 
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logarithmic factor. For a circular lamina of radius 1 and thickness 2b 
we take the demagnetizing coefficient v as that of an oblate spheroid with 
the same major axis and the same volume. The minor axis a of this 
spheroid is then 4b/7, and v=7a/4l=b/l. We will also require explicit 
expressions for s and A, which are : 


s=4l:; A=rla=A4lb. 
The expression for # than becomes 
L=—(1)n)[h/\/—b/l) -—Hev/(1— 4/6)... (9) 
Minimizing with respect to 6b we have : 


hbjl He Alb 
(be ~ (1— AJ) 


Assuming that /> b> 4, and expanding by the binomial theorem, we have 
b= (Aln) {1— sul V/(m)-3/ Vm) + O*)§ - se (1D) 


where p=(4/1)}/?. 
Substituting (11) in (9) and simplifying, we have 


L=—(1/n)H2[1 —n—pv/(q)— (3/9 —7)/8 +003) 2. (12) 
or to the first order 


L=—(He/n){(1—n)?—2n(1—)-V/(m)}- es (18) 


(10) 


§3. THE MAGNETIZATION CURVE. 


We shall now show that our theory gives a curve basically similar to that 
of Andrew (1948 b), based on the Landau (1943) theory. Initially the 
magnetization increases linearly with the applied field, and according to 
the zero order approximation discussed by Peierls (1936) and London 
(1937), this continues until 7=1—mn, whereupon the body enters the inter- 
mediate state and the magnetization falls linearly with further increase of 
yn, reaching the value zero when »=1. We shall see that, just as 
in Landau’s theory, the pure superconducting state actually persists for a 
short range above »=1—n to a value p, say, where the intermediate 
state becomes thermodynamically favourable ; the magnetization then 
falls abruptly to the value appropriate to the intermediate state. The 
condition for the transition is equality of the Landau functions for the 
superconducting and intermediate states. Since for the pure super- 
conducting phase 4 = H??/(1—n), this gives 


p*/(1—n)=1—(1/n){(1—p)?—2p(1 —p)V/(p)+O(H?)}. . (14) 


The zero order result p=1—n follows at once, and the next two approxi- 
mations are 


p=(1—n) {1+ pt?ny/(2)/(1—n)*4} 
and p=(1L—n) {1 yl ?2ny/(2)/(1—n)1/4—p3/4n21/4(2—n) 121, (15) 


Model of the Intermediate State of Superconductors 967 


The magnetization in the intermediate state then falls smoothly with 
increase of field. The magnetic moment / of the specimen is 


Gai 


~ i | Specimen < 4 Pall. S—domains a 
(for in the N-regions F= # ; in the S-regions F=0). 
Thus —4rM=hVpA/{(1—v)(1—npA)}, 


where V is the volume of the specimen. 


Hence —temS00[ (HSA) ], 
n l—y [l= 


or, after simplification, 


— 47 M=(HoV/n)[1—y—3u(39—1)/-V/(n) +O(2)]..- (16) 
It appears from (16) that the magnetization falls to zero where 
1—n=$pu(3n—1)/4/(m) or, since »<1, y=1—p. . . . . (17) 


But we observe from (11) that when 7-1, b remains finite and tends to 
(Al). This means that M cannot approach zero continuously, but 
has to do so stepwise ; for it is » which tends to zero, and the number of 
domains cannot be fractional. Thus ™ will fall to zero when the last 
S-domain is expelled. When there is just one domain 


—4r M = | H dr =sh|(1—v) (al? . 2b) 


one S—domain 
or, since / is nearly H,, 
—47r7M=H,. 2nl?b. 


Using this value in (16), and putting 7=1 in the terms in p, we have 
H2nPb=—4nM =(VHon){1——p+ OW2)} 
Whence 7=1—p—2z/?bn/V or, since b=1/(Al/n) or b/l=p, for 7 nearly 1, 
Gala bean). som. (18) 


Also the field for no domain is y7=1—p, by (17). Hence a linear inter- 
polation gives the field at which the last domain is expelled as 


Np=1—p(l +a). 5 ee ee Gare oe ERD) 


If the field is reduced, some hysteresis effects should appear, even under 
ideal conditions. For even apart from the “‘supercooling”’ of the normal 
phase (see the region P in fig. 3) there is the effect arising from the “ horn ” 
(region CBG in fig. 3). As there is no mechanism available for expelling 
the flux at C (fig. 3), the intermediate state will continue metastably to G, 
the place where (7) gives pA=1. This will occur for 


n= (1—n)/[(1— 4s/A)(1—v)] 
=(1—n)[1—$u(2—n)/1/(lL—n)+O(u?)]. . . - « (20) 
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$4. COMPARISON WITH EXPERIMENT. 
There are three sets of experiments to be discussed, viz., 
(a) Désirant and Shoenberg (1948) ; Andrew (1948 a). 
(b) Andrew and Lock (1950). 
(c) Meshkovsky and Shalnikov (1947 a, b) ; Meshkovsky (1949). 
(a) Only the experiments of Désirant and Shoenberg, who found the 
magnetization curves for long cylinders in transverse fields, need be 


considered, as Andrew’s resistance measurements under similar conditions 
give values of p in fair agreement with those of Déesirant and Shoenberg. 


Fig. 3. 


0 20 40 60 90) OER 
h (gauss)—> 


Magnetization curve for cylinder in transverse field. © increasing fields; 
+ decreasing fields. The experimental points shown are those for Désirant 
and Shoenberg’s tin specimen 89 at 3-:000° K. The theoretical curve is 
drawn for 4/r=0-67 x 10-3, i. e. d4=0-5 x 10-5 em. 


There are several ways of estimating 4 from the observations, but as it 
is (M/l)1/*, or in some cases even (4/l)'/4, that appears in the formule 
the magnetization curve is rather insensitive to the assumed value of A. 
Also the model is a considerable over-simplification, and a more rigorous 
theory may well involve extra numerical factors of order unity, which 
will be squared or raised to the fourth power if one attempts to calculate a 
value of 4 from an experimental result. Thus we may find, for instance 
that the value of » chosen to fit one aspect of the curve differs by a factor 
of about 3 from the best value for fitting some other aspect, and this would 
lead to a discrepancy of order 3? or 10 in 4. These difficulties also appear 
when one attempts to fit the branched model to the experiments. 
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Four different methods of fitting a value of 4 to the magnetization curve 
will be discussed. These are based on 


(i.) the value of p (fig. 3, BC). 
(ii.) the value of 7 where M has half the maximum value it attains for 
the bulk curve. 
(ili.) the value of 7 at M=0. 
(iv.) the slope of the falling part CD of the curve (fig. 3). 


We shall see that values of 4 obtained from (i.), (ii.) and (iii.) agree quite 
well with each other, but not with those from (iv.), which are too large by a 
factor of 50. 


Fig. 4. 


0-55 


Tile 10 


-4 
ale 20 30 40 x10 


Size variation of p for tin cylinders. © 3-5° K., assuming 4—0-55 x 10~° em. ; 
A 3-0°K., assuming A=0-28x10--cm.; [] 21°K., assuming 
A=0-19 x 10 em. 

(i.) Fig. 4 shows the theoretical curve for p as a function of 4/1. Owing 
to the slow convergence of (15), this curve has been obtained by solving 
(14) for » as a function of p. The experimental points of Désirant and 
Shoenberg are indicated, assuming the same temperature-dependence of 4 
as they did, but reducing their actual values of 4 by a factor 3-6. This 
makes 4,=0:55x10-5cm. at 3-5°K., 0-28x10-5cm. at 3-0°K., and 
0:19 X10-5 em. at 2:1°K. The present curve is seen to fit the experimental 
points rather better than that predicted by Landau’s branching theory. 
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(ii.) For the bulk curve M=4H,/47 for y=3/4 ; according to (16) this 
would occur in the actual specimen for 
$=2(1—y)—w(3q—-1)/-V(m) + OH"), 


since n=} for a long cylinder in a transverse field. Inserting the zero 
approximation =? in the second term on the R.H.S., we find the next 
approximation, 


or p= /(3)(8—4y)/5. 
Also from the geometry of the magnetization curve (fig. 3) 


A—Npe ee 3 


no—npn Me 2n¢ 
or 1—Np=(ne—np)/4ng= —(1/e—1)/4. 
Thus p= v/(3)(4—1/o—4np)/5, . ta ee 


where o=7q/np. 
The values of » and of 4 deduced from Désirant and Shoenberg’s experi- 
mental results by means of (21) are shown in the Table under the headings 
p, and 4;,. 

(iii.) According to (19), 7)»=1—2-6y, for a cylinder. The values of u 
and A estimated from Désirant and Shoenberg’s values of np are shown in 
the Table under the headings y;;;, and 4. 


(iv.) The slope of the magnetization curve at 7=3/4 is given by 


Heidi ae abe eee esa 


2{—1+23y/64/(3)+O(u?)} . 
From the geometry of fig. 3 


47, aM 2ng i Vv - rn 
He V dyn iad nD—Na (since —My=2n, . Ho V/4z7) 


=20/(1—c) , 
whence o/(1=o)=1--4(23/64/3) = "ne ee 


The values of » and 4 deduced from the experimental values of o by (22) 
are shown as p,, and 4, in the Table. 

The high values of 4 given by method (iv.) are perhaps due to the 
sensitiveness of this method to irreversible effects. For if (as in ‘fig. 3) 
the experimental values of M, for n just above p, are above the theoretical 
curve, the slope is likely to be greatly altered, while method (ii.) will 


give almost unaltered values of 4. The values from (ii.) are thus likely 
to be the most reliable. 
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Fig. 3 is a typical magnetization curve. Actually it is the curve for 
u—0-026. The points are those found by Désirant and Shoenberg for the 
specimen S9 at 3-000° K. For S9 /=75 x 10-4 em., and we have taken 
A=0:5x10-5. (This is a typical value, from 4;)). The curve is in 
slightly worse general agreement with the experimental points than the 
curve predicted by Andrew (1948 b) from Landau’s branched model. 
The experiments show a smooth fall of M to zero, rather than the stepwise 
termination we have suggested. The absence of steps may perhaps be 
explicable by using a thread model in the region where the ratio of super- 
conducting to normal volumes is small. eae 

(6) The experiments of Andrew and Lock (1950) on the magnetization 
curves of thin plates may also be compared with the theory in several 
ways. The consistency proves to be much poorer than in (a), and the 
general trend is to give very much larger values for 4. Little weight 
should, however, be given to the numerical values of 4, for the reasons 
given on p. 968 and also because we have assumed a laminar structure 
in which the lamine are circular, each of radius / (which need not be 


Flat plate showing one ‘‘ toroidal ’’ S-domain. 


constant for all Jamine). While this is a reasonable assumption for the 
cases of the sphere and the cylinder, the flat disk will probably form 
“ toroidal”? lamin (see fig. 5), and if we unfold them into flat lamin 
they will have the form of long thin rectangles, rather than circles. 
There is evidence for such toroidal domains in that the magnetization 
curves for falling field show large amounts of “‘ trapped” flux, suggesting 
the presence of multiply connected S-domains. It is very difficult to 
estimate the correction required to take account of this effect, but we 
might expect the value of 4 calculated on the basis of flat lamine to be 
too great by a factor of order \/(d/l), where d is the mean diameter of 
such a lamina, 7. e. of the order of the radius of the disk. This factor, 
of order 10 to 10?, might bring the results into reasonable accord with 
those deduced from the experiments of Désirant and Shoenberg. How- 
ever, comparisons which depend on the functional form of our expressions 
rather than on the actual value of 4 are less likely to depend on such 
details of the laminar structure and should give a more reliable check. 
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(i.) As in (a), we compare the values of p with the theory. The experi- 
mental work does not give the values of p directly, but rather J,,,., the 
maximum intensity of magnetization. As the initial slope “Si the 
magnetization curve is 1/(1—n) 


—4rI/h=1/(1—n) 
and in particular 
atlas opin)... 2 27 2 4. (28) 
Also from (15) 
p/(1—n)=1+-n4/(2)(1—n)-1/4 1/2, 


Fig. 6. 


~7 -—6 -5 -4 ~3 =2 
log 22y 


Logarithmic graph of 47Imq,/Hc—1 against 2lv. [-] 2-5°K.; A, 2:99°K; 
© 3-49° K. Full curve is theoretical. 


In these experiments n is nearly unity and, therefore, | 
Ar I paxfHo=1+-Ve2tAtiaa—ns. 2... (24) 


Andrew and Lock have plotted log {(—47Ing2/Hc¢)—1} against 
log 21(1—n) (21=L in their notation), and if the present theory is correct, 
the curve (fig. 6) should be a straight line of slope —4, in marked contrast 
to the theory of Landau and Andrew, for which the slope ought to be 
—2. It is seen that our law fits the experimental points very much 
better than the —? power law. It deviates significantly only for the 
thinnest specimens, where we might expect our approximations to 
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break down; in particular it is no longer likely to be true that we can 
take the lengths of the domains to be 2/ (since the regions outside the 
body from which the field is excluded become important, see fig. 2). 
The value of 4 deduced from this curve (fig. 6) is, however, very high: 
thus at 2/(1—n)=10-3 

anlage lls 
4./€. © “4h, nfl 2 V2) «Ae eee 


giving 4=3-5 x 10-* cm. 


Fig. 7. 


—3:5 -—30 log 9 (22) -25 —2-0 —!°5 


Logarithmic graph of o/(l—n)—1 against 21. A 2-99°K., © 3-49° K. 


(ii.) and (iii.) have not been attempted for these experiments since 
insufficient data exist. In particular, Andrew and Lock did not record 
the value of H;, in all of their experiments. 

(iv.) Here we have 


ise d (In M -1 5) 
rei taeereell —a—n)| 


from the geometry of the magnetization curve, where o has now been 
computed from the slope of the curve at 7=3. Hence, from (16), 


of(i— nel 4/(2)A a eee 


A noteworthy difference appears between the present theory and the 
branching theory, for where the latter requires o/(1—n)—1 proportional 
to (2/)-*°, our theory requires proportionality to (21)-1/. In figs 7 
log {a/(l—n)—1} is plotted against log (27) for these experiments, 
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and the slopes of the lines drawn are found to be —0:49 and —0-52 
at 3-5 and 3-0° K. respectively. However, the value of 4 deduced 
from the curve is again rather high: thus where 2/—10-3, we find 
log {o/(1—n)—1}=1-883 for the 3-5°K. curve. Hence 4Q=9-9 x 10-2 cm. 
at 3-5° K. 

(c) Shalnikov (1945) and Meshkovsky and Shalnikov (1948 a) proved 
the inhomogeneity of the intermediate state directly by observing the 
field distribution along a diameter in the equatorial plane of a split sphere. 
Later experiments (Meshkovsky and Shalnikov 1948 b) show that the 
non-uniformity of the field is also found on the surface of the sphere. 
The most recent work in this seriés of experiments (Meshkovsky 1949) 
summarized in the form of a contour map of the field in the slit between 
their hemispheres, shows the structure to be far from simple, both thread- 
like and laminar domains being in evidence. 

The size of the S-domains is directly measurable, using a case such as 
(A) or (B), fig. 11, p. 9, of their 1947a paper (7=3-16° K., h=72 Oe.), 
where the fraction of superconductor is small, so that our model of 
S-domains in a matrix of N-metal is as realistic as possible. In their 
fig. 11 (A), the number of domains is about 15, over a length of 4 cm., 
suggesting a domain size b=0-03 cm. (since x,, the fraction of N-metal 
=0-81); while a direct measurement of a typical single domain gave 
b=0-015 cm.; thus b=0-02 cm. may be regarded as a typical value. 
Substituting 4=10-°cm. (a rather high estimate) in equation (11), 
with J=4cm., we find b=7 x 10-cm., which is too small by a factor of 3. 
However, Meshkovsky (1949) comments on the apparent presence of 
a “fine structure” in the domains which would reduce the observed 
value of b, thus improving the agreement with the theoretical predictions. 

It is of interest to compare the values of b predicted by different models. 
Landau (1937) found b= ,/(7Al/2) for 7 near to unity for an unbranched 
model. This is similar to our form. On the other hand, for a branched 
model he found b=8 A?!*/1/3, which with 4=3 x 10-5 cm. gives 6=0-01 cm. 
Unfortunately no experiments to test the variation of domain size with 
specimen size have yet been reported. 

Meshkovsky and Shalnikov (1947 a) have reported a difference between 
the structures observed (A) as h/ is raised at constant temperature 7’, 
and (B) as 7 is lowered at constant h. The difference is probably not 
due to the fact that h is varied in (A) and 7 in (B), but simply to the 
fact that in (A) a,, the fraction of N-metal, is decreased, while in (B) it 
is increased; for two disjoint S-regions will exert very little influence 
on each other, and if they arise from two different nuclei (as is the case 
for a, increasing), there will be less correlation in their positions than if 
they arise from a single lamina which has broken up. 


I wish to express my thanks to Mr. H. N. V. Temperley, my supervisor, 
and to Dr. D. Shoenberg and other members of the Mond Laboratory for 
many helpful discussions. I am indebted to Cambridge University for 
a Research Maintenance Grant, during the tenure of which this work was 
done. 
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APPENDIX. 
By H. N. V. TEMPERLEY. 


Landau’s proof (1943) that L=fd7[F—(#B+hB—-h#)/87] is the 
appropriate potential to use in the theory when the field h is specified 
is incomplete. In particular he has not given a proof that L is the 
only function satisfying the required conditions. The potential, ¥ per 
unit volume, is assumed by Landau to satisfy 


OW ae. 
(a), 79 (Fr) » EB iat Aeeona 
together with h=nB+(1—n)H. ovena Patee 


Let Y=F-+ W, where F is the Helmholtz free energy per unit volume ; 
dF=H dB/4n—SdT. Then 


Megas 
q oF _ (oF oH = tee 
ae (Fr) (Se)n aH )ng=® C7) 
Therefore —(nB—h)/47(1—n)+ (s5) ==()s 
TB 
ow 
(sr) =(H—B)/4r7=(h—B)/4r(1—n), 
T,B 
and on integration 
W = (nB?/2—hB-+-h?/2)/4r(1—n). 
From (2) 
1/(1—n)=(H—B)|(h—B) ;_ n/(1—n) =—(h—H) [(h—B). 
Whence W=(—HB—hB+hH)/87, 
and Y=F—(HB+hB—-hH)/8r ..... ., (28) 


apart from an additive constant. 

As we require L={ Yd for the whole body, we may replace H, B 
by #, ZB, since the volume integration will average these variables, 
provided that the dimensions of the body are large compared with the 
domain size. We therefore have finally 


L=|[F-(HBAhB-hH)[8rj dr. we. (29) 
For a body in the intermediate state L takes the form 


Year+ |p. dS, ies (80) 


interfaces 


L= | Bs de | 


/ S-regions N-regions 
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where Ys; and Wy are the Landau functions per unit volume in the 
S- and N-regions respectively and y is the surface energy per unit N—S 
surface, introduced on p. 963. If Fy and Fy are the corresponding 
Helmholtz free energies per unit volume, then 


and Py=Fy—H/87 (since in N-regions =H.) Also if F%, is the 
Helmholtz free energy for the N-metal in zero field, 


» Py=FS AH? |8a 


= PF o+(H?2,4+-/?)/ 87. 
‘Thus 


Fodr+ | a eee | h(B—H) dr|8n 


N-+S-regions /¥ N-regions / S-regions 


ae | y.d8. 
interfaces 


=const-+ | H?, dz/8n—h (B—H) dt|8a+ | y ds. 
. » N-regions S-regions ~ interfaces 
This result involves no approximations in the penetration law. If, 
however, the S-domains are thick compared with the penetration layer, 
we may, as Ginsburg (1945) has shown, replace y by y—AH?./87= 4H?2,/87 
and omit explicit reference to the field penetration. Then 


L=const.+ | FH, do[80-+ | il dr] 8+ | (AH? /87) dS. 
N-regions S-regions interfaces 
(1) 
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SUMMARY. 


A calculation has been made of the consequences of assuming that: 
stars are produced in a cascade process by charged and uncharged particles. 
according to a very simple law. The constants in this law can be chosen 
to give agreement with experiment to the extent to which the quantities. 
are known and also fall within the ranges that are a priori permissible. 
It does not appear possible to determine them with much accuracy from 
existing data, but it seems that a fair proportion of the stars with less. 
than five prongs are due to primaries of fairly high energy. 


§1. INTRODUCTION. 


THE opinion has been growing for some time that the particles responsible. 
for the production of the stars seen in photographic plates exposed to 
cosmic rays can and do make more than one star before coming to rest- 
(Thomson 1949, Harding 1949, 195la). Also it is known (Camerini 
et al. 1949) that, in the larger stars at least, secondaries of large energy are 
possible which can themselves make stars. This opens up the possibility 
that the production of stars is a “‘ cascade ”’ process, and this view has been. 
developed in some detail by Messel and Ritson (1950) following a sug- 
gestion in a note to a paper by Heitler and Janossy (1949). Messel and 
Ritson did not, however, distinguish between the production of stars by 
charged and by neutral particlest. This has been studied by the Bristol 
school (Brown (et al. 1949, Camerini et al. 1949). They find that the- 
proportion of charged primaries is greater for the larger stars and is. 
greater high up in the atmosphere. It seemed of interest to see if the 
cascade theory can account for this effect and what, if anything, can be- 
deduced as to the production and energy of the secondaries. 

In order to do so we made the simplest assumptions which seemed to 
offer any chance of representing the main facts, since the mathematics 
of cascades readily becomes complicated and makes it difficult to see: 
the physics behind the calculations. 


* Communicated by the Authors. 
+ Messel has since done so (1951 and private communication) 
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Leaving aside the question of obliquity, to which we return in §7, we 
suppose particles moving vertically downwards through the atmosphere. 
These are of two classes—charged and uncharged. The latter will probably 
mostly be neutrons since the lifetime of the neutral mesons is too small 
for them to play any large part, and photons are usually excluded as an 
important factor because little or no correlation has so far been observed 
between stars and the “ showers ”’ of the soft component. Since photons 
of sufficient energy can cause stars in the laboratory, this conclusion: 
may prove too sweeping, but we shall accept it here. The charged 
particles will be mostly protons and z-mesons except right at the top of 
the atmosphere, where the true primaries still survive in appreciable 
numbers and «-particles, not to say heavier nuclei, become important. 
Because of this difference, our calculations will not apply to the top of 
the atmosphere. We shall assume that protons and 7-mesons have the 
same cross-section for star production (Camerini et al. 1950) and the same 
chance of producing secondaries. Most of the 7-mesons will decay before 
they make a nuclear collision; this simply means that the number effective 
for our purpose is less than the total actually produced. We take the 
cross-section for star-production, both by neutral and charged particles, 
to be the same, and denote the corresponding free path in gm./cm.? by X. 
In each act of star production we suppose that ¢ charged and ¢’ neutral 
particles, capable themselves of producing stars, are formed, each of which 
goes forward with the fraction qg of the energy E of the primary. Thus 
E[1—¢+¢’)q] is the energy used in those particles of the stars which 
have too little energy to produce further stars, and in producing z-mesons. 
which decay in the air. The lower limit for star production is written Ko. 
We assume an energy spectrum of the form N(E) dEa E™ to be valid 
at some point of the atmosphere which we take as origin, x=0. Though 
the form of our solution is such that this law will not hold strictly elsewhere, 
our results show (see §8) that the divergence is not great, so that it is 
not a matter of critical importance where the origin is taken. 

We suppose that the only way in which the neutrals lose energy is by 
making stars, but that the positives lose energy by ionization at the rate 
p per g./cm.? of path corresponding to minimum ionization. 

In §2 we give the mathematical solution of the cascade problem, 
in §3 we discuss the most suitable values of the quantities A, p, Dra wpe 
d, ¢ and n, in §§4, 5 and 6 numerical values are given, in §§7 and 8 
the effects of obliquity and of variations of the initial constants are 
considered, and in §§9 and 10 we set out the comparison between our 
theory and the experiments. 


§2. CascADE EQUATION. 


Let N (E, z)dE be the number of charged particles/cm.?/sec. at depth 
a gm./em.2 in the energy range E to E+dE capable of making stars, and 
let N’(E, x)dE be the same for uncharged particles. 


a x2 


980 Sir G. Thomson and P. E. Hodgson on the 
Then the rates of aie of these numbers of particles are given by 
ON(E, x) _ pe E ,(E =) | +2 an ON(E, «) 


and : 
ON'(E,x)_ ——N'(E z) | [s : e) LN’ (= “| | 
pee rt ae f 


chy A Aq q 
(1) 


The first term on the right-hand side of these equations represents the 
particles lost by star production and the second the particles produced in 
more energetic stars of energy E/g. The last term in the first equation 
represents the reduction in the number of charged particles due to 
ionization losses. 

We solve these equations by assuming a solution in the form of a power 
series : 


N(E, x) =exp (—a/A)E~" (a) +a, (2/A) + 9(2/A)?L. aes eee (2) 
N’(E, x)=exp (—2/A)E~"{bp +b, (%/A)+b9(2/A)?+. .. .}, 
where @,,a,....6,,b,.... depend on E but not on 2, and ay and by are 
constants. The exp (—a/A) factor is suggested by the exponential 
decrease in the number of particles with depth and the E~” factor by 
the power form of the energy spectrum. 

By substituting the equations (2) in the equations (1) one obtains, 
after some reduction, the relations 
24,0  3Agh* 


n—-1 
pL = iy Dy Pd, Be ee 
r r X ie 


d d 
toa it ag GG) +---}, ae oter(i) ae 


by 2byay  Bb,a" 
“ ¥ Se seal {0+ bo (+01) oe | 
BE E/q 
(3) 


where the suffixes outside the brackets show the energy at which the 
a’s and the b’s are to be evaluated. For convenience of calculation we 
write E=pdA/a, v=dq""1, v'=¢’q""1. 

By equating coefficients in the equations (3) recurrence relations may 
be derived : 


a 


1Q,=Ay_ V—Ap_..NX— oe et Vs te 
oe (4) 
rb,=v'(a yy +0n-4 ), 


where @,, 6, mean that a,,b, are reckoned at E/q instead of E. These 
relations enable the a’s and b’s to be calculated as power series in «. 
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§ 3. CHOICE oF VALUES OF THE CONSTANTS OF THE THEORY. 


In this section the most suitable values of the constants A, p, Ep; Y, 9; 
¢’ and n are discussed. | 

If the cross-section for star production is geometrical, then the inter- 
action mean free path A will be about 60 gm./em.?. But Bernardini et al. 
(1950) found a value of only 45 per cent of geometrical for the production 
of stars by 350-400 MeV. protons. For shower particles, which are 
mostly z-mesons, Camerini e¢ al. (1950) found a cross-section of 95-39 
per cent of geometrical for production of secondary stars in emulsion. 
Although Bernardini’s results refer to much lower energies than we are 
considering, we choose a value of A nearer to his than to Camerini’s for 
two reasons. First, we are mainly concerned with the effects of protons 
as he was, and not of z-mesons, which produced most of Camerini’s stars ; 
secondly, we propose to consider only stars of five and more tracks, since 
the smaller stars can be produced by particles of less energy then we can 
consider, but some of the smaller stars are likely to be produced by high 
energy particles, making grazing impacts, and hence the cross-section 
for the larger stars can hardly have the full geometrical value. We 
take a value. for the cross-section of one-half the geometrical, which 
corresponds to A ~ 120 gm./cm.?. 

The rate of energy loss for a relativistic proton is about 1-5 MeV./gm./ 
em.2 (Braddick 1939). The lowest proton energy of which account is 
taken in this theory is 450 MeV. (see below), and these have a rate of 
energy loss only 30 per cent above the minimum. It is therefore a good 
approximation to take the ionization energy loss of the charged particles 
as constant. 

The series (2) converge only if «7/A<1, and since the atmosphere, when. 
obliquity is taken into account, is about 12 mean free paths, | 2/A| may 
reach the value 6 and we must have «<%. In practice, for reasonably 
rapid convergence, we take ace. Since pa Ph TOK" this means 
that E>1800 MeV. Particles of greater than the limiting energy 
E,=1800 MeV. are referred to as of high energy. 

In addition, it is possible to calculate the number of particles in the 

‘energy range from E,qg to E, MeV. (see §5). Particles in this range are 
referred to as of medium energy. Particles of energy less than E,q, 
the slow particles, are not covered by this theory. 

The value of g, the fraction of the energy of the primary given to the 
average secondary, is taken as ¢. 

The lowest energy of particles considered here, E,qg—Kp, is, therefore 
450 MeV., which corresponds on the average to stars with five or more 
heavy tracks (Barton, George and Jason 1951). 

The stars produced by high energy particles are mainly of more than 
ten tracks, for which Brown et al. (1949) found about one shower particle 
per star at 11,000 ft., while Camerini et al. (1949) found 2-0 per star at 
70,000 ft. Many of the stars of ten or more tracks but no shower particles 
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will in fact have been produced by medium energy particles, and, if all 
these stars are omitted, the corresponding numbers of shower particles 
per star at the two altitudes are 2-3 and 3-2. Also, Fowler has shown 
that 75 per cent of the shower particles are 7-mesons, which mostly decay 
in flight before producing a star. These particles must be excluded by 
dividing the above figures by four. The value of ¢ is consequently rather 
uncertain, and we take ¢=%. This is probably a slight over-estimate. 

Harding (1949, 1951a) has shown that in most cases the primary 
particle when neutral retains a large fraction of its energy. It will count 
as a secondary for our purposes, and there may of course be others. 
Allowing these to compensate the cases where the primary does not 
retain enough energy to make another star, we take ¢’=1. 

Probably ¢ and ¢’ vary rather considerably with energy, whereas here 
we take them independent of it. This is unfortunate, but for the higher 
energy particles the theory only depends on ¢ and ¢’ through v and v’. 
Now while ¢ and ¢’ probably increase with energy, ¢ falls. Thus v and 
v* vary more slowly with E than ¢ and ¢’. 

The exponent n in the energy spectrum, both for primary particles 
(Hoang 1951) and for stars (Barton et al. 1951), has been shown to be about 
2-5, and we took this value for convenience of calculation, though some 
recent experiments indicate that it may be a little high. 


§ 4. Higu ENerRGy ParTIcuEs. 


From (2) the number of high energy charged particles at depth 2 is 


fi exp (—a/A)E~?° {a-ha (3) +... 4 dE 


+E, 


exp (—2/A) p21 zs 
Semin inee si (4045+ co faze 1) da. 


The coefficients a,,@).... b,,b,.... were calculated as power series in 
« multiplied by a, or by and the integral evaluated for various 2/A. 
A similar calculation was made for the high energy neutral particles. 


§5. Meprum Enrerey Parricuss. 


The medium energy particles at any depth x are produced in nuclear 
interactions above x. It is therefore necessary to evaluate the number 
produced at any height z that reach x and integrate from the top of the 
atmosphere down to x. 

The rate of production of neutral particles of energy E to E-+dE 
produced at a height z is 

¢ dE 

7 LN(E/g, 2) +-N'(B/¢, lore . 
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The fraction that reaches a depth « is exp [—(x—z)/A]. Therefore the 
total number of medium energy uncharged particles at x is 


ve Dee 
[|S oxp -@—2y a} NGBla, 2)+N'(H/g, 2] aE ae 
-HJ x, AY 
where H is the height of the top of the atmosphere from the origin of x. 
‘This was taken at —7A. 

This expression can be integrated exactly, and was evaluated as a 
function ofz/A. The corresponding calculation for medium energy charged 
particles has to take into account the energy p(w—z) lost as the particle 
goes from z to a The expression for the total number of medium 
energy charged particles is 

rx rE, +p(4—2) 


kes ie Ha yee {—(@—2)/A}IN(B/q, 2) +N(E/q, 2)] dE dz. 


This cannot be integrated exactly, and so graphical methods were used. 


§ 6. RESULTS OF CALCULATION. 


In this section the results of the calculations of the last two sections 
are presented. These results for each class of particles consist of a part 
in a) added to a part in by), where a, and by are adjustable constants. 
‘The relative numbers of the various types of particles depend on the 
ratio do/bp. 

The variation of the neutral and charged particles and of the high 
energy and medium energy particles with atmospheric depth 2/A is 
plotted in figs. 1 and 2 in arbitrary units for various d/b). The factor 
-exp (—a/A) is omitted from the curves shown as its inclusion would make 
it harder to see the run of the curves. 

The variation of the ratio of charged to uncharged particles with 
atmospheric depth is plotted in fig. 3 for various values of d/o. 


§ 7. THE ANGULAR DISTRIBUTION OF THE STAR-PRODUCING RADIATION. 


The theory as developed so far is one-dimensional, and takes no account 
of the angular spread of the cosmic radiation. 

If it is assumed that the secondaries all continue in the same direction 
as their primary, then it is possible in principle, though rather laborious 
in practice, to extend the theory to three dimensions. Instead, we take 
the angular spread of the cosmic rays into account by altering the scale 
of the curves by a factor equal to the average ratio of the vertical distance 
between two atmospheric depths and the distance covered by a particle 
in going from one depth to the other. This is only an approximate 
‘method since the scale factor alters with atmospheric depth. Its use is 
equivalent to assuming that the angular distribution of the cosmic 
zadiation does not vary with depth. 
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The radiation at the top of the earth’s atmosphere is roughly isotropic, 
but, as it proceeds through the atmosphere, the particles whose directions 
make a large angle with the vertical are preferentially absorbed, with the 
result that the radiation becomes more collimated as it penetrates the 
atmosphere. But the secondaries do not continue in the direction of 
their primaries, having an angular distribution about them. This tends. 
to make the radiation less collimated. 


Fig. 1. 
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Variation of intensities of charged and neutral particles with depth. 


These two effects oppose one another, with the result that over the middle 
and lower regions of the atmosphere the angular distribution does not 
vary much. This conclusion is supported by the work of the Bristol 
group (Brown et al. 1949) and of Lord and Schein, whose results on the 
angular distribution of the star-producing radiation at 11,000 ft. and at. 
53,000 ft. show little difference. 

These results give an average obliquity factor of 0-8. The assumption 
of constancy of angular distribution, while probably approximately 
correct for the middle and lower regions of the atmosphere, is certainly 
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quite incorrect near the top. The results of this theory are not therefore 
applicable to the first one or two collision mean free paths of the atmo- 
sphere. This would have been the case even if the complication due to 
obliquity was absent, owing to the presence of the heavy primaries. 

This method of allowing for obliquity is undoubtedly crude, but a 
more refined treatment would have to take into account not only the 
angular distribution of the star-producing radiation but also the angular 


Fig. 2. 
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Variation of intensities of medium and high-energy particles with depth. 


distribution of the secondaries about the direction of motion of the: 
primaries, which introduces great mathematical complexity. Itis preferable 
to using the Gross transformation and neglecting the latter effect, since 
this predicts a considerable increase in collimation. 


§ 8. Errect oF VARIATIONS OF THE INITIAL CONSTANTS. 

0 / 
The constants whose values are the most uncertain are ¢, ¢’, q and A. 
If there were no cascade or obliquity effect, the absorption mean free 
path as observed should equal the interaction mean free path, which we: 
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have assumed to be 120gm./cm.2.. The effect of the cascade is to increase 
the expected absorption mean free path to 150 gm./em.°. Correction for 
obliquity brings it back to about 120 gm./em.?, which is close to the 
‘observed value. 

The absorption mean free path of 150 gm./cm.? is thus composed of a 
part depending on A and a part depending on ¢, ¢’ and q. The same 
-absorption mean free path can be obtained for a smaller value of A by 
increasing the latter constants. But a limit is set to this process by the 


Fig. 3. 


Ratio of charged to uncharged particles 


Variation with depth of the ratio of charged to uncharged particles. 


‘approximate experimental values of ¢ and ¢’, and by the inequality 
(P+¢')q<1, where the inequality must allow both for the slow particles 
and the stars and most of the 7-mesons. It is not therefore possible to 
vary very much the value of \ we have taken. 

Increasing both ¢ and ¢’ increases the ratio of charged to uncharged 
particles, since it makes the relative effect of the ionization loss less 
important. The ratio of charged to uncharged particles is also increased 
by increasing ¢/¢’. 

When ¢ and ¢’ have been chosen, the above inequality and the known 
‘characteristics of stars sets limits on the value of g. The main features 
of the theory do not depend critically on the choice of q. 
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§ 9. VARIATION OF THE RATE OF STAR PRODUCTION AND OF THE 
ENERGY DISTRIBUTION OF THE STAR-PRODUCING 
RADIATION WITH HEIGHT. 


The variation of the rate of star production with height is shown in 
“fig. 4 for various ao/by. It is closely exponential, giving absorption 
“mean free paths of 145, 150 and 155 gm./cm.? for a /b)=1, 4 and 4 


Fig. 4. 
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Rate of star production in arbitrary units 


120 gm fem? 
116 gm/cm? 


= Gir 4 ee 2h a0 2 4 6 


Variation of the rate of star production with depth. 


‘respectively. The corresponding mean free paths when the angular 
-distribution of the star-producing radiation is taken into account are 116, 
120 and 124 gm./cm.?. 
The energy distribution of the star-producing radiation was calculated 
-at three heights, and it was found to be still closely a power spectrum, 
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but with slightly different index. The index of the power law increases 
with altitude from 2-2 at sea level to 3-0 three-quarters of the way up 
the atmosphere. 


§ 10. CoMPARISON OF THE RESULTS OF THE THEORY WITH EXPERIMENT. 


The calculated ratio of charged to uncharged particles decreases with - 
atmospheric depth, in accord with experiment, for Camerini et al. (1949), 
found a ratio of 0-58 for stars of five or more tracks and 0-47 for stars 
of four or more tracks at 70,000 ft., while Brown et al. (1949) found. 
corresponding ratios of 0-36 and 0-26 at 11,000 ft. (The slight minimum 
in the calculated curve is not considered significant.) The value at the 
higher altitude can be fitted by suitable choice of a/b), but that at the 
lower altitude is higher than predicted by the theory. However, Brown 
et al. (1949) find only 10-8 stars per c.c. per day at the Jungfrau, while 
Harding (1951 b) finds 16-2. Assuming, as is likely, that the missing stars 
are mostly small, for which the ratio is very low, the values 0-36 and 0-26 
will be too high. 

In making these comparisons it should be remembered that the theory 
refers to a cascade in air, whereas the stars investigated by the Bristol 
group occurred in photographic emulsion. In the latter medium, unlike 
air, charged z-mesons frequently produce stars, thus increasing the 
apparent ratio of stars with charged to stars with uncharged primaries. 
This effect cannot be evaluated, since the 7-meson flux depends on the 
amount of matter surrounding the plates, which is not known, but taking 
it into account would improve the agreement between theory and 
experiment, 

The calculated ratio of medium to high energy star-producing particles. 
at 11,000 ft. is 6-1, and on our scheme only the high energy particles will 
be able to produce shower particles. Since d=, and only one in four 
of the shower particles are protons, this gives 0-37 shower particles per 
star, in good agreement with the values 0-38 and 0-29 found by Brown 
et al. for stars of five or more and four or more tracks respectively. 

Our calculated exponential variation of the rate of star production has. 
been found by many investigators (see Messel and Ritson 1950). The 
observed absorption length varies with energy and, therefore, with star 
size (Addario and Tamburino 1949) and falls from about 145 gm./em.2: 
for star sizes 3, 4, 5 to about 113 gm./cm.? for star sizes greater than 6, 
which agrees sufficiently with the calculated range of 116-124 gm./em.2. 

Laboratory observations show that particles of a given energy do not all 
produce stars of the same size, but give quite a broad distribution of 
star sizes. This makes a comparison between theory and experiment 
difficult. 

Our theory does not take account of the small stars. These may be 
due either to particles of low energy or to particles of any energy which 
make glancing collisions with the nucleus. For the present it seems 
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difficult to disentangle the two classes, but the small value we are led 
to assume for the eross section for star production is an indication that the 
second class is important. 
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XCIX. The Deflection of a Cable due to a Single Point Load. 


By E. Markianp, 
The University of Nottingham ™*. 


[Received May 28, 1951.] 


SUMMARY. 


A method is given for calculating the deflected shape of a uniform: 
inextensible cable supported with each end at the same level, when a 
small point load is placed anywhere along the cable. 

The analytical treatment is straightforward and a method of successive- 
approximations for the numerical solution of the equations is suggested.. 
Results have been computed for the case of span to sag ratio of 10 and- 
ratio of point load to cable weight of 0-1, and are compared with. 
experimental values. 


THe method of flexibility coefficients has recently been introduced’ 
(Pugsley 1949 a and b) for analysing a stiffened suspension bridge. To- 
obtain a numerical solution it is necessary to know the behaviour of a 
suspended cable when a single load is placed anywhere along it. Professor~ 
Pugsley emphasized that the deflection of any point of the cable is not a 
linear function of the applied load, and the required deflections for an 
arbitrary magnitude of applied load were obtained experimentally. 

The calculation of the deflected shape of the cable, although presenting" 
no theoretical difficulty, is arithmetically troublesome. It has been found, 
however, that the problem may be solved reasonably quickly by the- 
method given below. 

Consider the uniform, heavy, flexible cable of weight W shown in 
fig. 1, in which the ends A and B are supported at the same level. It is. 
assumed that the initial dip at the centre, D,, is small compared with the- 
span, 2L, so that the initial shape ACB is very nearly parabolic and the- 
load per unit horizontal length is constant. 

When a vertical load w is applied to the point Go, the cable deflects. 
to the shape AGB, AG and BG being parabolic arcs. 

The horizontal distance of the loaded point from the mid-point of AB. 
increases from aL to aL and the vertical distance below AB from D, to D.. 
The horizontal movement is taken to be small so that the loads per unit 
horizontal length of the ares AG and BG are equal. 

2HL 


Let SA ese : 
e als aM ee a na ek OAS ce THN 


* Communicated by the Author. 
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where H is the horizontal component of cable tension. Then the equations 
of AG and BG referred to the rectangular axes shown on the figure are 


: 2 
for AG :— y=Awt 5, Wao eade Nteeeayie stows (2) 


PA 
for BG :— y=Aget >, Men Oe ees iste et (3) 


where A, and A, are the slopes of the arcs AG and BG at G. Resolving 
forces vertically at G, we obtain 
2L w 


A,+A,= rcanw.. . . . . . . ° (4) 


Substituting x=—L(1—qa) in equation (2) and ~=L(1+a) in equation (3) 
and solving for A, and A,, we obtain 


ib 2 
for AG :— Me = | a+9 F +a | a+ = rae ret ry Pe 13) 
ii 2 
for BG :— y= | a sy —«| e+ = pe | een O)y 


For a specified value of a), the lengths of arcs AG, and BG, of a given: 
cable may be computed. Let these be 8, and 8, respectively. Then, 
since the cable is inextensible, 8;= are length AG, S,=are length BG. 


Hence 
iguteal L2y2\1/2  ¢ Lu ( wey u=(L+aq-+1 
es epee me sana 1 
L 2 u(t C ) +f bee} C « a C } u=(L+ay +a 
(7): 
and 
w 
Soect T2a2\ 1/2 ¢ Lu : ae u=(l-a)y +1 
Ts =F u(14 ce + 7 log C + + C2 epee 


‘992 E. Markland on the 


By inserting the known initial values of S,/L, and 8,/L into these equations, 
we have two simultaneous equations in a and c/L which may be solved 
by trial for a given value of w/W. The numerical work is helped by 
forming approximate values for 8, and 8, as follows : 


L(1—a) ihe ay, 19 
S.= | {14+ 5[ (Hayy +a)+2 | ae de, 


‘L+a) 2 fa gy a a2 


Upon integration, and subsequent partial differentiation with respect to 
.a and c/L respectively, the following results are obtained : 


oS Lies 27 w /—w 2 wr {Ss Ne 
tia los w( Faye +1) —(14+ wy) (Tage ta) |, 
— __ O(8,/L) _ L?2 [w (——w 2 w\ (—— w 2 
tart Cote Ae oat be (may, +1) (1+ 9) ean -a) ] 
ao OSL) geet CLNe w 3 (—w 3 


peated 10) (mate +ip-(mag—)']. 


Fig. 2. 


O/%2L 


O:24%D- > 


ve 


Approximate values of these partial differential coefficients may be 

‘computed by inserting trial values of a and c/L, then the total differentials 
d(S,/L)=«, da+B, d(c/L), 
d(S_/L)=a, da-+By d(c/L) 

are known approximately. 

If trial values of a and c/L yield discrepancies 4(S,/L) and A(S,/L) 
im equations (7) and (8), the required corrections 4a and A(c/L) are 
-given approximately by 

ja B24SuL)—B, AS 9/L) 
%1B,—o 2B, 
and (10) 
A(c/L)= eS (aco eae) : 
%1B,— oP, 
By making corrections in this way it is found that equations (7) and (8) 
may be solved after a few attempts. 
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The equations of AG and BG are therefore known and the difference 
in ordinates between them and the original curve AG,CB may be computed. 
Strictly, the vertical component of displacement of a point on the cable 
is not the difference in ordinates since each point of the cable has, in 
general, a horizontal displacement as well. However, tests on a cable 
in which horizontal deflections were measured at each of nine points in 
turn, show that the maximum horizontal deflection at an unloaded point 
was 0-10 per cent of the span for 2L/D,=10 and w/W=0-1, and was 
measured at a point one-fifth of the span from one end. The maximum 
error in the calculated vertical deflection is, as indicated on fig. 2, about 
one-quarter of one per cent of the central dip. Both horizontal and 
vertical components of the displacement of the loaded point are of course 
obtained correctly from the calculations. 

Calculated results for the nine stations defined on fig. 3 are given below : 


Fig. 3. 


Station 


TABLE. 


Vertical deflections expressed as percentage of central dip. 
‘Span: dip ratio=10: 1. w/W=0:1. 


: Load at station no. 
Deflection at 


station no. l 2 3 4 5 
Ls +4-32 42-59 +0-82 —0-47 —1-35 
2 42-59 | “oa, +2-57 +0-09 —1-63 
3 +0-79 £5952 +4-99 +1-68 —0-85 
A —0-59 —0-01 +1-65 +4-24 4102 
5 SO Ray ae She —0-93 +0-97 +3-95 
SS bei, | yeas Zogo) 1-29 +1-02 
Lad eae —3-25 ere eae Hee a 
Goan A188 —2-89 —3-13 2272 «| 2168 
9 e195 —1-81 —2-03 —1-88 Sage 
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Tests were made on a chain of about 5 feet length and 0-93 lb. weight 
and deflections were measured at each of nine stations along the chain with 
a travelling microscope. Some trouble was experienced, especially 
towards the ends, with friction between the links, but the results shown 
on fig. 4 agree reasonably well with calculated values. Experimental 
values previously published appear to be somewhat at variance with the 
present results. 


Fig. 4. 


Load at station no. / 


u © A Load at station no. 3 
Load at station no. 5 


+ y The curves are pletted from the 
2 calculated valves of TableL 


DEFLEXION +INITIAL CENTRAL DIP % 


~3 


STATION NY 


n 


a 


° Load at station na 2 


+} x Load at station no.4 
3 

The curves are plotted from the 
py 


caleviated valves of Table { 


° 


4 


y 
N 


DEFLEXION + IN/TIAL CENTRAL D/P Yo 


-3 


STAT/ON NP 


As a check on the accuracy of the apparatus, the deflection of the 
central point has been measured over a wide range of loads applied at 
this point, and a comparsion with calculations in which the simplifying 
assumption of a parabolic initial shape is not made, is shown on fig. 5. 
The method of this straightforward calculation is given in the Appendix. 

The span to sag ratio of the experimental chain was varied between 
5 and 10, and it is interesting to note that deflections expressed as a 
fraction of the central dip varied little over this range. It appears, 
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therefore, that the values given in the Table would not vary greatly for the 


different ratios of span to dip likely to be encountered in suspension 
bridge practice. 


Acknowledgment is due to Mr. R. A. Dorton for his help with the 
experimental work and to Mr. R. C. Coates for many valuable suggestions. 
The author is indebted to Professor J. A. Pope for placing the facilities 


of the Engineering Department and Workshop of Nottingham University 
at his disposal. 


DEFLEXION AT CENFRE 
Load applied at centre 
© Experimental results 


The curve «s based on the 
method of Appenaix A 


DEFLEXION = INITIAL CENTRAL D/P % 
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Referring to the centrally loaded cable of fig. 6, the equations relating 
the coordinates of any point of the cable AG to the intrinsic angle at that 
point are : ; 
WiC==SOCAVNE cle nates an hes es 1 (LD 
aje=log (sec o+tan ys), . . . « - (2) 

3x¥2 
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and the distance along the curve from P is given by 


sfe=stan Wh ews) 4) <> en ee 
Hence 
Dje=sec ,—sec fg, . 2) so 6) ee 4) 
sec #,-+tan oy 
L/c=log (Se) o. 4s a eRe re (5) 
and 
S_ tan #,—tan pq (6) 
Li 5). /secsb s+ tan shy \ ee re 
> ( %gttan wb, 
Fig. 6, 


x Oo 
Now 
tan by w+W 
tan bq re eas 
or 
Ue tan bq (7) 
W. -tanwd,—tandag? 9 5 eae 


The computing procedure adopted was to assign a value to %, and to 
calculate %, from equation (6) above knowing the value of S/L. D/C 
and L/C follow from equations (4) and (5), The corresponding value of 
w/W is obtained from equation (7). 
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C. The State of Ionization in Lithium Hydride. 


By M. 8S. Aumep, 


University College, London*. 


[Received May 7, 1951.] 


SUMMARY. 


Single crystal photographs of LiH have shown that the structure is of 
the NaCl-type, a=4-09 A., Z=4. The intensities of reflection from planes 
having all even and all odd indices give atomic form factors for Lit* and 
H~*respectively which, when extrapolated to (sin@)/A=0, lead to the value 
X=0:25+0-25. This result is not modified by taking into consideration 
the thermal vibrations of the atoms. It follows therefore that in the solid 
state LiH has only 25 per cent ionic character, with limiting values 
(much less probable) of 50 per cent ionic character or complete covalency. 


§1. InTRODUCTION. 


In the NaCl-like structures, it is well known that the percentage of ionic 
character may vary according to the polarizability of the component ions. 

The state of ionization is particularly interesting in the case of lithium 
hydride, since it might at first be thought that when Li and H combine 
chemically to form LiH, a complete transference of an electron from the 
outermost L shell of Li to the innermost shell of H would occur; thus 
both Li+ and H~- ions would have the helium stable configuration. In 
this paper it is shown by X-ray methods that this is not the case for LiH 
in the solid crystallized state. 

The structure of lithium hydride has been determined by different 
workers, using the powder method. Many discrepancies appear in some 
of the work that has already been done. Some research workers have, 
by extrapolation of their results, concluded that LiH possesses 100 per 
cent ionic character. Other assume without proof, that in LiH both Li 
and H ions possess two electrons. 


§2. ReporTED WoRK ON THE STRUCTURE AND STATE OF IONIZATION 
in LrH. 


(A) Bijvoet and Karssen (1922) prepared powder photographs for Li 
using Cu and Cr radiations. Because of the chemical reactivity of the 
substance, photographs of LiO and LiOH were also taken in order to 
eliminate lines due to them arising from the possible decomposition of 


LiH. Their results were as follows : 
on ee BL I a RA a a pee BA ol eo da oa eee 
* Communicated by Professor Kathleen Lonsdale, F.R.S. 
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1. The observed lines were in the positions which would be expected 
from a NaCl or ZnS arrangements, with a unit cell edge of 410kX; four 
molecules per cell. 

2. They then made a comparison between observed and calculated 
intensities (see following Table). The first column gives the indices of 
the reflections, columns 2 and 3 the observed intensities using Cu and 
Cr radiation respectively, columns 4 and 5 the calculated intensities 
assuming a zincblende and NaCl arrangement respectively. If Li has 
the NaCl grouping, then, they argued, it would be anticipated that 
+-ve Li ions and —ve H ions would be present in its crystals. 


Indices T,ps(CuKa) Tops(CrKa) Teatc(Zn8-type) — Leare( NaCl-type) 
111 medium weak weak 10 0 
200 strong medium strong 0 10, 
220 medium strong medium strong 9-5 8-8 
113 medium very strong 6-5 0 
222 very weak strong 0 3°6 


From the above Table they concluded that there is a better fit with the 
NaCl arrangement, than with the zincblende arrangement. 

3. The poor agreement obtained they explained by saying that it is 
not to be expected that with such simple atoms there would be good 
agreement with calculations which neglected the electronic distribution. 

(B) Bijvoet (1923) published a further paper on the structure of Li and 
Lit but with no improvement in the observed intensities of powder lines 
for LiH. He mentions that after the exposures in most cases the lithium 
hydride had been reduced 15—20 per cent by weight. The time of exposure 
varied from 2 to 12 hours. In some of the long exposures the sample was 
renewed after a few hours. ‘ 

Now it is clear that the intensity values given cannot be reliable in 
these circumstances. The relative intensities, using Cu and Cr radiations, 
do not seem to be consistent among themselves or to agree particularly 
well either with the NaCl-like or with the zincblende-like arrangement. 

(C). In 1929, Bijvoet and Frederikse published a paper called ‘‘ The 
Scattering power for X-rays and the electron distribution of the H~ ion ” 
in which curves were given for (f,;+f4) and (f,;—f_) (but without any 
experimental intensity measurements to support them). From these 
curves it was argued that complete ionization is almost certain, but it 
was admitted that the photometrically measured values of the intensities 
(not given) showed considerable fluctuations amongst themselves. 
Observations of intensities of reflections appearing on single crystal 
rotation photographs do not agree with the data given by these workers. 
It is probable that decomposition of their LiH powder occurred during 
the exposures. 

(D). Zintl and Harder (1931) also determined the structure of LiH 
using the powder method. The unit cell constant was given as 
4:084 kX (+0:001). These workers obtained reproducible values for 


5 
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the intensities. CuK« radiation was used and the intensities were measured 


by a photometer. The intensity formula used, neglecting absorption and 
temperatures factors, is 


1+ cos? 26 
r Fre 
sin’ 6 cos 9° ? af 
CZ Dieu (p’=multiplicity). 


The intensities given by Zintl and Harder are :— 


hkl Tynot pz EZ co vs 
Lai 11-0 64-8 ‘17 
200 > 20-5 34:5 ‘60 
220 10-5, 29-5 36 
311 7-0: 41-3 Ath 
222 3:5 12-7 “28 
400 1-5 8-2 ‘18 
331 3-0 39:0 ‘O09 
420 4-0 37-1 ait 


From the above Table they concluded that (F5,_»,)?>(Fs,-9n 41)? 
always, and that therefore Lill has the NaCl-like structure. Ifthe structure 
were of the zincblende type the observed structure factors would have to 
satisfy the relations 


(Bean) Espagne)" span a) 


In their discussion on LiH, Zintl and Harder pointed out that if Lit 

and H~- ions existed, each with two electrons, then all the “odd” 
reflections should disappear. 
.. They claimed that the appearance of all the odd reflections can be 
explained as being the same effect found by Debye and Scherrer (1918) 
in NaF, by Gerlach and Pauli (1921) in MgO, and by Mark and Tolksdorf 
(1925) in SrCl, : namely a differential decrease of scattering power f of 
different sorts of atomic particle, with increasing 0. 

They pointed to a reported change in the conductivity of LiH. near the 
melting point, which although smaller than, is similar to that for alkali 
halides, and they concluded that the conductivity of LiH is cationic. 
‘This does not necessitate, however, that Li and H are completely ionized 
in the solid state of LiH. 

Assuming complete ionization in LiH, Zintl and Harder calculated the 
radius of the H~ ion to be equal to 1-26 kX units (using r,;+ 0-78 as 
given by Goldschmidt). 

(Z). In 1932, two short notes, one by Bijvoet and Karssen and a 
reply to it by Zint] and Harder were published : but these did not clarify 
‘the position further. Had Zintl and Harder drawn the graphs (f,;+/1) 
and (f;—f) for the “ all even ” and “all odd ”’ reflections respectively, 
thus deducing the graphs for f,; and fy, they could have seen that Li and 
H are by no means completely ionized, as they assumed them to be. 
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Since that time it seems to have been taken for granted that LiH is 
completely ionized. The following work will show that this cannot be 
the case. 


§3. SINGLE CrysTaL RoTaTION PHOTOGRAPHS. 


By careful cleavage from a rather big crystal of LiH it was possible 
to obtain small crystals without distortion on any of the crystal faces. 
These were quickly sealed into Lindemann glass tubes, and no unwanted 
powder lines were found on the short-exposure rotation photographs 
required for intensity measurements. Powder lines showing orientation 
were, however, observed on rotation photographs of larger crystals used 
for ‘‘ diffuse spot ’” measurements. 

The crystal was set with [110] vertical, and Cu radiation was used from 
Victor and Machlett tubes run under similar conditions. Intensity 
measurements were made using spots appearing on the zero layer line 
only. 

The crystal dimensions were roughly 0-25~0-30mm. along [100], 
0-35mm. along [010], and 0:35~0-40 mm. along [001]. The crystal, 
while rotating, was completely bathed in the incident filtered X-ray beam. 

Cylindrical rotation photographs were taken, with times of exposure of 
one, two, four, five and ten minutes, and the multiple film technique was 
used. 

With such short exposures the intensity of the incident X-ray beam 
was fairly constant, and white radiation radial streaks could not be 
observed on the rotation photographs, which showed small undistorted 
Bragg spots on a clear background. 

With the Weissenberg camera, an intensity scale was made for eye 
estimation of intensities, using the 004 reflection from the same crystal. 

The unit cell dimension from rotation photographs is a=4-086 A. 
The value given by Zintl and Harder (formerly expressed in kX units) is 
a=4-092 A. 


Observed structure factors. 


In the following Table, the observed mean relative intensities from 
different rotation photographs are given for comparison with those of 
previous observers. The observed structure factors are obtained after 
allowing for the multiplicity, polarization and Lorentz factors only ; 

329 
z denotes ct , p the multiplicity on the rotation photograph. 

The single crystal rotation photographs of LiH, therefore, confirm that 
the structure is of a NaCl type. 

In Table I. the observed structure factors F, obtained from the single 


crystal method, and those given by Zintl and Harder using the powder 


method, F,,, are compared. 


It may be seen that the ratio F UP, is nearly constant, within the 
experimental errors of both methods. 
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Atomic scattering factors. 


Atomic scattering curves for the Lit* ion and for the H~* ion (uncor- 
rected for temperature vibration) may be deduced from the observed 
structure factors curves for “all even’’ reflections (fiuitfa) and “ all 
odd ” reflections (f;;—fq). For any given value of (sin )/A 


Fut iutSa)+ (fu—fa)}, 
fat Vus-le)—(fu—fa)}- 
The extrapolation of the graphs for f,; and fy to sin 6/A=0 is difficult 
because of the wide gap involved, but taking the intercepts on the 
sin 6/A=0 ordinate as OY, OY’ for Li and H respectively, it is found that 
possible limits of the ratio OY/OY’ are 3: 1 and 2-5: 1-5. 
Now this intercept is proportional to the number of electrons in the 


atom. Hence OY/OY’=3—X/1+X with 0-5>X>0. The most probable - 
value is X=0-25. 


This means that the Li-H bond has 25 per cent ionic character. The 
limiting values would be 50 per cent ionic character or complete covalency. 


TABLE [. 
sin 6 Te i (sind)? 
hkl = —— p.2 Tpnot (ws oF?) F,=futfa F, F/F, eB(") Fs 
111 ee OD 6-7 2-6 41 6:3 -93 2-8 
002 245 2-16 60 27-8 5:3 77, 6:8 Ol 5:8 
ye 8 3465 1-3) 25 19-1 4-4 60 7:3 83 5:3 
ihe: 406 2-12 14 6-6 2-6 41 62 To 33033 
222 424 2-04 20 9-8 3:1 3) | a!) 765 4-] 
004 490 1:03 8 7:8 2°8 42 46-6 “69 4-() 
331 534 2-4 8 3°3 1-8 30 «61 -465 2-8 
224 600 42 15 3°6 1-9 ‘571 3:3 
115-333 -636 18-04 30 1-66 1-3 53D 2-4 
for any value ee » Fsp—2n i8 always > F sp -:en+1- 


§4. THERMAL VIBRATION IN LiH. 


The thermal scattering of LiH and LiF has been studied by means of 
“* diffuse ’’ Laue photographs and compared with that of KCl and other 
alkali halides (Ahmed: Thesis, London, 1950). The results suggest that 
the elastic behaviour of LiH is similar to that of LiF, but that the values 
of the elastic constants are higher for LiH than for LiF. The elastic 
behaviour of LiH is different from that of KCl. 

From the results of Giinther (1920), Ubbelohde (1936) calculated the 
Debye characteristic temperature for LiH and found 0)=815° K. at 
293° K. But @, is a measure of the restoring force on the displaced 
lattice points and so is larger the greater this force ; in other words, the 
harder the crystal. For simple crystals the characteristic temperature is 
a convenient measure of hardness in this sense (James and Brindley 1931). 
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0, for LiF is 648° K. at 273° K., which means that Li is harder than 
LiF. This is in accordance with the evidence of the diffuse spots them- 
selves that the elastic constants for LiH are greater than those for Lif’. 
The mean square amplitude of atomic vibrations in Lil at 293° K., eet 
for @, the value of 815° K., is given by Lonsdale (1948) as u?—0-0585 A?. 
The Debye-Waller formula, I=I,e-™, which applies to cubic crystals 
composed of one kind of atom only, at temperatures not near to the 


Fig. 1. 


x 0 


Atomic scattering curves uncorrected for thermal vibration, 


O1 02.03 O04 O05 06 O07 > 


melting point, gives F=F,e™ where F, is the structure factor for the 
ideal crystal at rest, F is the observed structure amplitude and 


2 
2M= its 


sin @\2 — sin@\ 2 : S722 aoe 
3 (>) . 2s YOR F=F,e3 (5) where B=—. u?. 


r 3 

In the case of LiH, a different value of B ought to be applied to each atom, 
since using the same value of B implies that Li and H should vibrate with 
equal amplitudes. The amplitude of vibration of the lighter H atom is 
probably greater than that of the Li atom, and in general, thermal 
vibrations must be different for every crystallographically different atom 
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in a unit cell, since they depend on the surroundings of the atom as well 
as on its inertia. Due to the uncertainty in the value of B however, to a 
first approximation B may be taken to be the same for Li and H in LiH. 


877? 


B= w=1-54 A? 


In Table I. the structure amplitudes F% corrected for the temperature 


factor e-® ey are tabulated in the last column. 

Thus the curves representing (f,;+/), and (f;;—/_) (corrected for the 
temperature factor) for the ‘“ all even ” and for the “ all odd ”’ reflections 
are obtained. From these, the separate atomic scattering curves for 
Li*~* ion and H-~* ion, (corrected for temperature scattering), are deduced, 


Fig. 2. 


27: 


V25b._ 


0 oO! O02 03 04 OS O68 O07 one 


Atomic scattering curves corrected for thermal vibration. 


and are shown in fig. 2. These give the same value of X as before. A 
relative variation in the temperature corrections for Li and H respectively 
would appear mainly in the high orders and would not affect the conclusions 
drawn concerning the state of ionization in LiH. 
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ABSTRACT. 


Numerous growth spirals with shapes ranging from circular to regular 
hexagonal have been observed on carborundum crystals of types I. and IL., 
by coating the crystal faces with a thin film of silver of reflectivity nearly 
90 per cent and then working in reflection, using phase contrast 
illumination. The shapes of these spirals, in relation to the crystal 
structure, are shown to be in accordance with the predictions of Frank’s 
theory. On a type II. crystal, hexagonal interlaced spirals have been 
observed. The annihilation of growth fronts where they meet, and other 
properties of growth fronts are illustrated and hence the observed patterns 
for two or more screw dislocations are explained. The observed density 
of dislocations ranges from small values up: to a maximum of 
104 per sq. cm. For a typical circular spiral the calculated values of the 
radius of the critical nucleus is 24 and supersaturation is 0-2 per cent. 
By the application of multiple-beam interference, using both Fizeau 
fringes and fringes of equal chromatic order, the step heights in spirals 
have been measured accurately for type II. crystal and found to be 
15 A. This is equal to the height of the unit cell, proving that these are 
growth spirals originating from screw dislocations, in exact accordance 
with theoretical prediction. 


§1. INTRODUCTION. 


Accorpine to Burton, Cabrera and Frank (1949) growth of crystals at 
low supersaturations can take place only in the presence of dislocations. 
A crystal should, therefore, have a number of dislocations with a screw 
component terminating on the crystal face. Frank (1949) showed that 
when growth takes place on the molecular terraces so exposed, the edges 
of monomolecular layers develop as “ growth spirals’ centred on the 
dislocation. Experimental evidence supporting this was found on beryl 
by Griffin (1950). The present paper deals with dislocations and “ growth 
spirals” found on the faces of carborundum crystals. 


* Communicated by Professor 8S. Tolansky. 
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§2. Tum CrystaL STRUCTURE AND TYPES OF CARBORUNDUM. 


Carborundum, a strongly homopolar crystal, occurs in at least eight 
known types. One of these is cubic, sometimes called B-Sic. All other 
types (a-SiC) are based on either hexagonal or rhombohedral unit cells. 
All types have identical layers, but differ in their arrangements. Each 
type is distinguished by the number of layers in the unit cell. The basal 
pinacoid is predominantly developed, and it is on this that these growth 
spirals have been observed. 

The crystals studied were either pale green, dark green or black. Some, 
which have a few well developed faces, have had their structural types 
determined by a goniometric method (Thibault 1944). X-ray diffraction 
methods will be required for the others. 

No cubic crystal has been studied. Most of the crystals examined were 
of the commonest type, carborundum II. Its unit cell is hexagonal 
with six layers, and in Ramsdell’s notation is 6 H (Ramsdell 1947). The 
lattice parameters are 


a=3-073 A, c=15-079 A. 


In this work the only other type available was carborundum I. This 
has a rhombohedral unit cell and referred to the hexagonal axes consists 
of fifteen layers. In Ramsdell’s notation it is 15 R with 


a=3-073 A, c=37-70 A. 


§3. VISIBILITY OF GROWTH SPIRALS. 


99 


Numerous “ growth spirals’? have been observed on the faces of 
carborundum (SiC) erystals. These spirals are well developed ; over 
thirty turns of the spiral can be traced in some cases. The spirals are 
centred on a dot which must mark the point of emergence of a screw 
dislocation on the crystal face. 

The clean surface of the crystal when examined by a metallurgical 
microscope, using bright field illumination, does not show up these 
features, the surface appearing smooth. A little amount of impurity 
makes these features slightly visible (fig. 8, Pl. XXVI. has been taken 
by this method). Breathing lightly on the crystal face, while it is under 
observation, increases the visibility and the “ lines ” flash out, but soon 
disappear when the water re-evaporates. The “‘ lines’ seen by breathing 
are dotted and slightly diffuse, but the visibility obtained by this method 
is high. Though convenient for visual observations, the procedure is not 
suitable for photography. 

Following Griffin’s technique, contrast was increased by deposition of 
a thin film of silver on the crystal face using thermal evaporation. This 
was improved by the use of phase contrast illumination. The 
photographs have been taken using positive phase contrast, the absorption 
of the phase plate being 80 per cent and phase retardation equal 
to $7. With this equipment many features are often still faint and only 
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just visible. Usual photographic methods for increasing the contrast 
have been employed in some of the photographs given. 

In fig. 13 (Pl. XXVIII.), it is seen that the visibility and contrast is high 
compared with other photographs. The surface of this crystal was not 
very clean when it was silvered. The silver deposit was not uniform, 
the surface appearing mottled. Examination with phase contrast 
illumination gave high visibility. This suggests that impurity can make 
the steps more visible. 


§4. SHAPE OF THE SPIRAL FOR A SINGLE Screw DisLocaTIon. 


The growth spiral will have a shape depending upon the rate of advance 
of a growth front in different crystallographic directions. 

For growth taking place in accordance with Frank’s ideas two cases 
arise :— 


(a) when the Frenkel kinks on the step (exchange sites) are close 
together and the distance moved by an adsorbed molecule before 
it hits a step is large ; 

(6) when the kinks are few or the distance moved by the adsorbed 
molecule is small. 


When the former conditions exist, which are more likely to occur in 
growth from vapour, the molecule will have a high probability of adhering 
to the step if adsorbed near it, irrespective of the crystallographic 
orientation of the step. Thus for growth from the vapour, taking the 
rate of advance of the ledge to be independent of direction, the ledge 
will form a simple spiral under steady uniform supersaturation. The 
form of the spiral can be represented by the Archimedian equation 


f= 20.0, 
and the constant spacing between turns is 
ér=4rp,, 


where p, is the radius of critical nucleus, equal to a¢/2kT In«, and where 
is the supersaturation ratio, a is the interatomic distance and ¢ is the 
neighbour-neighbour binding energy of the crystal. 

These predictions are confirmed by the circular spirals shown in 
figs. 1, 2, 3 and 8 (Pls. XXV., XXVI.). Itis to be noted, however, that the 
spiral turns are more closely spaced at the centre, and the spacing gradually 
increases on going away from the centre until it becomes nearly constant. 
This is in accordance with the prediction of theory (Burton, Cabrera 
and Frank 1951). 

Any dependence of the rate of advance of a growth front on the 
orientation of the step-line could impose a distortion of the growth 
spiruis, so that they exhibit the symmetry of the crystal face. Such an. 
eftect has been found in fig. 4 (PI. XXVI.), where the spiral is nearly circular 
at the centre, and gradually takes on the hexagonal symmetry of the 
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crystal face (carborundum type II.). In fig. 7 (PI. XXVI.) is a hexagonal 
spiral observed on a crystal identified to be of type II. This is a regular 
hexagonal spiral, showing clearly the dependence of the growth rate on 
the crystallographic orientation. The step-lines very near the centre 
show a curvature, which, as predicted by theory, decreases with distance 
from the centre. However, soon the straight step-lines start rounding 
off at the corners and become gradually curved on moving away from the 
centre. These step-lines do not meet the line of discontinuity sharply, 
but tend to curve away from it, as shown in fig. 7, for the outer five or 
six step-lines. 

The straight edges show that under certain conditions there is a sharp 
minimum in the growth rate as a function of orientation. | 

In figs. 10 and 11 (Pl. X XVII.) further spirals are shown. Fig. 10 has 
been observed on a rhombohedral crystal type I. ; 


$5. PROPERTIES AND BEHAVIOUR OF GROWTH FRONTS. 


The growth pattern for two or more screw dislocations ending on a face 
is complex and depends on the properties of the growth fronts. The 
growth fronts starting from a single screw dislocation will spread on the 
surface of the crystal face in regular shapes, but the presence of a 
boundary or an obstruction will distort its regular shape. When the 
advancing growth fronts meet an obstruction they can propagate round 
corners, as shown at the upper end of the obstruction in fig. 4; this 
point behaves almost as a source of secondary growth fronts. 

The growth fronts emitted by two sources of opposite sign annihilate 
each other where they meet. Numerous examples of this can be 
seen in figs. 2 and 3 (PI. XXV.). When there are several screw dislocations 
actively emitting growth fronts, a point on the crystal face will be in the 
dominant field of only one of these, and this alone determines the number 
of growth fronts passing through it. This is exemplified in fig. 2 where 
there are five screw dislocations. By choosing the point of observation, 
say, near the bottom screw dislocation, the number of growth fronts 
passing through it is solely determined as if this alone were active. 


$6. GROWTH SPIRALS FOR Two oR More Screw DiIsLocations. 


In fig. 17 is drawn schematically the growth pattern for two screw 
dislocations A and B of the same hand. The ledge starting from A goes 
on rotating and generating the spiral until it meets at p, the ledge 
originating from B where it terminates. At p, the two ledges fuse with 
one another. As the spirals unfold themselves the ledges meet at 
Pv Ps Ps » - . in the upper half and at ps, py, p, . . . in the lower half 
of the figure. The locus of points of intersection for two equal spirals 
has been shown to be a cartesian oval (Burton, Cabrera and Frank1951). 
Between p, and py (ps, 4), ete., the missing parts of the two spirals are 
shown by the dotted lines. The resultant figure is the solid line curve, 
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Starting from A and going round the spiral we descend by one step in 
each complete turn till we come to the point p,. From p, we go round the 
spiral B to the point p, and then again on the spiral A to the 
point p;. Thus in going round the resultant curve once, we descend by 
two steps from p, to p; to py and so on. Similarly, starting from B and 
going round we will descend by two steps, from p, to p,, etc... . As 
the curves gradually smooth out at the points of contact, the figure will 
appear to be two spirals alternately spaced. Fig. 17 has been drawn for 
' circular spirals, and the case for spirals with straight edges can similarly 
be drawn where the behaviour will be the same. In fig. 3 (Pl. XXV.) 
the case for circular spirals is illustrated. 


Fig. 17. 


Fig. 8 (Pl. X XVI.) illustrates a unique example in which the spiral is 
doubled with the members strictly ‘“‘ parallel’. It may be that the twin 
form arises from two close-by dislocations of the same sign, cooperating 
closely. 

In fig. 9 (Pl. X XVII.) is illustrated the case of two screw dislocations of 
the same sign close to each other. The spirals are circular and developed 
for two or three turns after which the ledges from other neighbouring 
dislocations interact with them. ice 

In fig. 10 (Pl. XXVII.) is shown another case of two similar screw 
dislocations close to each other, so that growth starts from both of them. 
The spirals are circular at the centre for one or two turns after which they 
settle down as spirals with straight edges. The meeting of the growth fronts 
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from the two screw dislocations results in small kinks in the innermost 
straight edges. As in the case of circular spirals, here also by going 
round we descend by two steps, and the figure has the appearance of 
two spirals alternately spaced. However, the step height between 
successive lines in any particular direction is just one unit. 

Fig. 11 (Pl. XX VII.) shows the case in which the central screw dislocation 
has developed several turns before reaching the point where the second 
small screw dislocation emerges. They quickly accommodate each other 
and the behaviour is similar to the previous cases. The central screw 
dislocation is dominant and would determine the rate of crystal growth 
of the face. Near the edge of the figure it is seen that some more screw 
dislocations emerge on the crystal face, and are situated in the same 
direction from the central screw dislocation as the second one. These 
new screw dislocations yield within a short while to the domination of 
the central one and quickly fall in line, passing on the growth fronts with 
slight delay and in slightly modified form. 

Now consider two screw dislocations of opposite hand terminating on 
the crystal face. The ledge starting from one terminates on the other, 
and with the unfolding of the spirals the two arms join together generating 
closed loops. This is exemplified in fig. 3 (Pl. XXV.), where one of the 
spirals has developed several turns before the other starts growing. 

The growth pattern for three screw dislocations of the same hand 
observed on a crystal of type IT is illustrated in fig. 12 (Pl. XX VIT.). Here, 
again, the spirals are circular at the centre and soon change into hexagonal 
spirals. It can easily be seen that, for three screw dislocations 
cooperating with each other, by following any one of the edges, the 
resultant curve descends three units in each turn. It should again be | 
emphasized that the step height in the successive lines will be just one unit. 

At the centre the ledges are closely spaced and, once they have settled 
down, the ledges are more widely spaced. 

Fig. 2 (Pl. XXV.) illustrates the interaction of five screw dislocations. 
The growth fronts emitted by the screw dislocation at the top of the 
figure meet, first, the growth fronts advancing from the screw dislocation 
on the left of the figure and then, later, the growth fronts advancing 
from the bottom screw dislocation. Annihilation occurs over portions of 
the resultant growth fronts which assume characteristic shapes composed 
of ares of circles with these dislocation points as centre. Near the centre 
of the figure, due to the interaction of different growth fronts a crowding 
occurs. Here the compound growth fronts are not ares of circles with 
the different dislocation points as centres, but instead there is a gradual 
change in curvature. 


$7. HexaGonaL INTERLACED SPIRAL. 
These features observed on a crystal face of carborundum type II. are 
shown in figs. 13 and 15 (Pl. XXVIII.). In each of these hexagonal spirals 
the edges fork out at the corners and meet the two neighbouring edges. 
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This interlacing gives the hexagonal spiral the “ spider web ” structure 
shown clearly in fig. 13 which shows the interlinking of four such 
Spirals. 

A double interlaced spiral is shown in fig. 14 (Pl. XXVIII.) in which the 
edges do not meet each other. The lines are broad with inner edges sharp 
and the outer diffuse. 

An explanation of this interlacing has been given by Frank (1951). 


$8. MEASUREMENT oF Step HEIGHT. 


The step height has been measured by the use of multiple-beam 
interference fringes using both Fizeau fringes of equal thickness and the 
fringes of equal chromatic order (Tolansky 1948). 

A rounded spiral, effectively a spiral conical hill, may be considered 
as a circular conical hill made of a series of concentric steps with a constant 
spacing d and step height h. It can be shown that if a Fizeau fringe 
passing over the peak of such a hill appears to bend through an angle 20, 
and if X is the dispersion, 7. e. the distance between successive fringes, 
then the step height h is given by 


h=dx sin 6/2X. 


‘By measuring X, 6 and d, h can be evaluated. Fig. 5 (Pl. X XVI.) shows 
the Fizeau fringes passing over the spiral of fig. 4 (A5461). 

For examination with fringes of equal chromatic order, the image of 
the peak of the hill is projected on to the spectrograph slit. This was 
arranged to prodce a system of fringes of equal chromatic order parallel 
to the spectrum line in one half of the field and sloping towards the red 
in the other (fig. 6, Pl. XXVI_.). 

Since the step height is obtained from the mean of all those contributing 
to form the peak, irregularities of the optical flat and the lack in 
faithfulness of silver contouring are averaged out. 

In one case another determination of the step height has been possible. 
In fig. 4 (Pl. XX VI.) starting from the centre of the spiral we can reach any 
point A on the line of discontinuity, either by going round the 
“ fish-like ” obstruction in a clockwise or an anti-clockwise direction. 
In the former, we encounter, say, N, steps. Therefore the point A is 
N, steps lower with respect to the centre of the spiral. In the latter case 
we go down, say, N, steps from the peak of the hill in order to arrive 
at A. Thus at A the difference in level on the two sides of the line of 
discontinuity is that corresponding to (N,—N,) step heights. This 
difference is quite large (over 50) and is therefore capable of being 
measured accurately and easily either by Fizeau or fringes of equal 
chromatic order. Knowing this, the step height has been evaluated, 
but is subject to uncertainty due to the possible presence of hidden 
dislocations in the obstacle and discontinuity surfaces. 


3822 
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$9. Resutt oF MEASUREMENT OF STEP HEIGHT. 
It is seen below that the step height measured in the cases so far is 
equal to just one unit cell. 


on EEE eee ee 
i 


Fringes of Method by - 
Fizeau fringes | equal chromatic | counting No. Crystals 
order of steps 
152A 771A 14-8 A ‘Spiral in 
fig. 4, Type II. 
141A 
145A Circular spiral, 
15-5 A right of fig. 2. 
15-0 A 
15:3 A 


§10. Discuss1IoN OF RESULTS. 

From fig .2 (Pl. X XV), which is at magnification x 90, the spacing between 
the successive arms of the circular spiral when it is nearly constant is 
approximately 2-5 mm. Using 6r=47p, we have 

(oy == 90 x4n min As, 2p. 
From p,=a¢/2kT In « and using Trouton’s rule to estimate ¢/kT according 


to which 4¢/kT=3-5Tb/T=6 at an absolute temperature of 0-6 of the 
boiling point of the material we get 


a —~, 0-2 per cent. 


$11. Denstry oF DisLocaTIons. 


The density of dislocations varies widely on different specimens of 
carborundum. In fig. 4 (Pl. X XVI.) there is only one which dominates 
the growth of the face. In figs. 1, 2, 3, 11 and 12 (Pls. XXV., X XVII.) 
there are several. The largest density of dislocations observed is shown in 
fig. 16 (Pl. XXVIII.). Here there are nearly 104 screw dislocations per 
square centimetre. 

Another noteworthy point is that in fig. 16 (Pl. XXVIII.) there are 
nearly twenty screw dislocations, most of which are of the right-handed 
type. It is characteristic that in any one region there is a large 
predominance either of right-handed or of left-handed screw dislocations. 


§12. Summary. 

Numerous growth spirals starting from points of smergence of screw 
dislocations have been photographed on the surfaces of carborundum 
crystals. These spirals have characteristic shapes. The shape is discussed 
in relation to the crystal structure. The behaviour and properties of 
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growth fronts and the growth patterns for two and more screw 
dislocations are illustrated. From the observations on the photographs 
it is shown that the type of information that can be derived is the 
determination of the radius of the critical nucleus, the supersaturation 
and density of dislocations. 

_ The step height has been measured accurately by multiple-beam 
interference fringes. The step height at the edges of these growth spirals 
is found to be one unit cell. 
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ABSTRACT. 


A discussion is given of certain features of the growth of carborundum, 
as illustrated in Verma’s paper. It is shown that polytypism may occur 
as a result of growth through the rotation of a spiral growth step. 


$1. INTRODUCTION. 

In the preceding paper, A. R. Verma describes growth markings that he 
has observed on carborundum, and shows that they are in accord with 
the theory of crystal growth developed by Burton, Cabrera and Frank 
(1949, 1951). This note deals with certain features, mainly peculiar to 
carborundum, which were not predicted in the paper quoted. I would 
like to thank Mr. Verma for showing me his photographs before 
publication. 


§2. SHAPE OF SPIRALS. 


According to our theory, growth of a crystal from the vapour should 
usually produce round growth spirals rather than polygonal ones. This 
conclusion follows from the supposition that similar forces bind a building 
unit into the crystal and into an adsorption position on a close-packed 
face: whence calculations show that on the closest-packed edge of a 
crystal monolayer the distance between “ Frenkel kinks” is likely to 
_ be small compared with the average distance which an adsorbed 
molecule diffuses on the surface. The edge can only advance by the 
adhesion of building units to these kinks; thus the rate of advance of 
a step is independent of its orientation and round spirals result. If, 
on the other hand, the distance between kinks is large, the spirals become 
polygonal, tending to be straight in the orientations of slow advance, 
which are those in which this distance is greatest. 

Carborundum, however, does not grow from its vapour. The most 
likely presumption is that the transport of C and Si from carbon and 
silica to the growing carborundum crystal proceeds by gas molecules of 
CO and SiO ; the oxygen being removed as CO, (or perhaps O,) which 
regenerates CO on the hot carbon, while CO also reduces SiO, to the 
relatively volatile SiO. Most probably the surface of the SiC crystal 
has a layer of combined oxygen, during growth, on which adsorbed 
SiO and CO molecules undergo surface diffusion. The state of the 
surface would then vary with temperature and with the CO-CO, balance 


* Communicated by the Author. 
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in the atmosphere. Thus the surface conditions would be much more 
complex than those we assumed. The appearance of both round and 
polygonal spirals, while still unexplained in detail, does not, therefore, 
constitute a failure of theoretical prediction. There is indication, 
e.g. in fig. 12, that the rounded form of spiral belongs to the final stage 
of growth, at which it might be connected either with a lower temperature 
or a less strongly reducing atmosphere. 


§3. VARIABILITY OF SPACING. 


In nearly every case the turns of the spirals are more closely spaced 
near the centre than further out. This must be attributed to increased 
supersaturation shortly before growth ceased, which would be the 
natural accompaniment of a fall in temperature, but which could again 
be complicated by changes in composition of the atmosphere. The 
variations are in some cases gradual, in others, e. g. figs. 11 and 12, abrupt. 
In fig. 11 the result is that previously dominated dislocations have begun 
to develop their own little growth pyramids of two or three spiral turns 
in extent. Before the increase of supersaturation the growth steps 
attached to these would never have developed a full turn before the 
arrival of another step from the dominating centre. 

Spacing is also liable to be non-uniform across a crystal face, as well 
as radially on each hill, as, for example, in fig. 11, where the widest 
spacing, indicating the fastest advance of steps, occurs in the right hand 
central region. The most probable cause is screening of the crystal 
surface by other crystals, though internal strain could also be responsible. 
The same cause must explain the dominance of the central dislocation 
over similar ones below it. 


$4. Ho~ttow DIsLocaTIons. 

In most of Verma’s pictures there is a plainly visible dot at the site 
of the dislocation, suggesting that its core is a hollow tube with a 
diameter of the order of magnitude 1 or 2 microns. This is generally 
to be expected from theory (Frank 1951) with dislocations of Burgers 
vector exceeding about 10 A, because the large strain energy near the 
dislocation is then insufficient to overcome the surface energy, which 
tends to close the tube up. 


§5. CRYSTAL STRUCTURES OF CARBORUNDUM : 
STACKING OPERATORS AND ZHDANOV SYMBOLS. 

Before proceeding further with the discussion of Verma’s photographs, 
it is necessary to consider the crystal structures of carborundum. In 
each of the known structures there is a silicon lattice, and an identical 
interpenetrating carbon lattice, which is one of the infinity of structures 
resulting from the layered close-packing of spheres. (The word layered 
is inserted here, because no mathematical proof appears yet to have been 
given for the strong presumption that every close-packing of spheres 
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can be dissected into close-packed layers.) As is well known, in all such 
structures the layers projected normally on to a plane parallel to one 
of them fall into three possible positions, by which they may be labelled 
a, b or c. The possible close-packings correspond to any sequence of 
these letters in which no two successive letters are alike. Taking any 
such sequence, for example, abcacb... , and letting these letters 
represent Si layers and the corresponding Greek letters represent 
C layers (or conversely), we obtain from it a carborundum structure 
ax bB cy axcy bB ... This places each atom of one kind at the centre 
of a tetrahedron formed by four of the other kind. 

We obtain a simpler notation by introducing a pair of “ stacking 
operators’, as suggested by Nabarro (private communication). We 
shall use 4 and V, where he proposed + and —. These describe the 
relationship between a pair of layers. . 6 after a, c after b, a after c (or, 
in carborundum, bf after ax, etc.) is denoted by 4: c after a, 6 after c or 
a after b is denoted by V. With suitable choice of coordinates, 4 or V 
' represents the orientation of a triangle of spheres providing a crevice 
into which a sphere of the next layer fits: and the notation has mnemonic 
value for the fact that a rotation of 60° converts the operation 4 into V. 

In the special case of carborundum, 4 and V have an alternative 
interpretation. Atoms of the crystal have strong lateral bonding in 
double layers («b), (Bc), etc. Instead of calling 4 and V stacking operators, 
in this case they may be considered to denote such layers directly, 
A signifying («b), (8c) or (ya) and V signifying (yb), (Ba) or (ac). It would 
not be a certain inference, however, that the surface layers during growth 
are necessarily of («b), etc., rather than (ax), ete., type : 7. e. whether the 
outer layer of atoms belonging to the crystal are triply or singly bonded 
to it. In the presence of gases containing oxygen this remains an open 
question. 

In this notation the structure ax bB cy ax cy bB becomes 44 4 VV. If 
the latter sequence is continuously repeated we have the structure of 
carborundum I. In abe notation the sequence only repeats after 15 
layers. The notation can be further contracted in Zhdanov’s symbol 
(Zhdanov 1945), which consists of pairs of numbers in each of which the 
first is a number of consecutive ds and the second the following number 
of consecutive Vs. We _ shall write it Z(DiPo PsPy - - > Poq-1 P2q)- 
Carborundum I. corresponds to Z(32). The commonest carborundum is 
carborundum II., Z(33), but there are many rare forms such as 
Ramsdell’s 87 R, Z(33 33 33 33 32). Cubic carborundum, Z(10), which is 
also rare, is the only one whose conditions of formation appear to be 
distinctive with any certainty. It has long been recognized that the 
relationship between the various forms of carborundum is so close as to 
deserve a special name “ polytypism ” (Baumhauer 1915, Thibault 1943). 

Ramsdell (1947) made an important generalization concerning the 
types of carborundum, from the eight types of ascertained structure 
at that date, Z(10), Z(22), Z(32), Z(33), Z(34), Z(33 32), Z(33 33 32), 
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Z(33 33 33 33 32): namely, that the number pair 33, which applies to 
the commonest type of carborundum, also occurs within the Zhdanov 
symbol of most of the others. He also noted the non-occurrence of 
the number 1. (Since Z(10) can also be written Z(nO) it does not make a 
necessary exception.) 

It seems probable that the Z(33) structure is one either of real 
thermodynamic stability, or of particularly favourable growth rate, and 
certain that some other cause than thermodynamic stability must be 
responsible for structures differing from this by a change in the packing 
at large regular intervals. This is clear enough with the last two structures 
listed above, without needing to refer to such extreme cases as one which 
repeats its layer sequence at intervals of 198 layers (Honjo et al. 1950). 

All the layered sphere-packings, and hence all the carborundum types 
have trigonal symmetry by simple rotation. One layer has 2-, 3- and 
6-fold axes normal to its plane. When another layer, or more, is placed 
on it, the 2-fold axes fail to coincide and vanish, while 6-fold axes 
coincide with 3-fold axes, so that 3-fold axes remain. If the Zhdanov 
symbol is divisible into two similar parts, each of an odd number of terms, 
e.g. Z(11), Z(33) or Z(31 1311), the structure becomes hexagonal by 
virtue of a 6, screw-axis. The crystal then has macroscopic hexagonal 
symmetry ; but the force field of any close-packed layer still only has 
trigonal symmetry. 

Apart from these hexagonal structures, the single cubic structure, and 
aperiodic structures, all remaining structures can be described as 
rhombohedral. When the difference between the number of 4 and 
V operators is zero or divisible by 3, however, there is an alternative 
trigonal cell containing one-third as many atoms: there is no known 
example in carborundum*. The rhombohedral unit cell extends through 
three times as many layers as the true stacking sequence: so that the 
structure Z(33 33 33 33 32) which repeats its pattern after 29 layers is 
called 87 R by Ramsdell. 


.§6, DisLocaTIONS AND POLYTYPISM. 


Polytypism as observed in carborundum is readily comprehensible 
in terms of the dislocation theory of crystal growth. A crystal only 
increases in thickness (save in the large supersaturations generally absent 
after nucleation of a sufficient number of crystals) by reason of the 
dislocations it contains. The crystal that grows hereafter is not composed 
of an indefinite number of layers stacked upon each other, as ideally 
considered, but of a finite number of interleaved helicoids (or other, more 
complex, expanded Riemann surfaces). Each helicoid axis is a dislocation. 
The structure necessarily repeats with a period corresponding to the 
pitch of the screw. 
pee ee ee ee ee ee 

* Note added in proof.—One such, Z(32 23), has now been discovered : 
Ramsdell, L. S. and Kohn, J. A. (1951), Acta. Cryst., 4, 111. 
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When one observes under the microscope the growth from aqueous 
solution of CdI, or of PbI, to which some Cdl, has been added (Forty 
1951) one may see very thin crystal plates increasing to a considerable 
lateral extension (say, | mm.) with perfectly uniform and constant 
interference tint, i. e. not increasing in thickness at all. Quite suddenly 
a group of growth pyramids will appear, followed perhaps by a second 
and third group, each appearing suddenly. After the growth of these 
pyramids has proceeded for a while the observer can begin to resolve 
spiral growth terraces on these pyramids. Ultimately one or two will 
dominate the growth and the crystal will continue to increase in thickness 
in a steady way by the slow spreading of turns of these spirals. It is 
thought that the initial plate-like growth, for high degrees of super- 
saturation, may occur by surface nucleation, without dislocations : and 
that the crystal then probably becomes self-stressed through the 
non-uniform distribution of impurities, ultimately up to its theoretical 
yield stress, when the thin plate buckles and shears, so raising terminated 
steps in its surface. Similar phenomena may well occur in initial thin 
tabular crystals of carborundum, though in this case there will be 
alternative sources of stress, for example, due to crystals in contact, and 
thermal stresses due to partially screened intense radiation, etc. When 
this happens, if the shear is by a uniform amount with a fairly abrupt 
termination, the structure which grows thereafter will be that of the 
exposed edge at the step. If the original crystal has some simple initial 
structure, say, Z(33), and slip is by an integral multiple of the repetition 
period of this structure, the same structure then grows in thickness. 
But if slip is by a non-integral multiple of this period the result is a 
polytype whose period equals the step height. Such slip leaves a misfit 
on the slipped surface, so should be less common than the integral 
slip. This would account for Z((33),, 32) structures. It can also generate 
the cubic structure. -It seems improbable that this always arises 
secondarily in this way from hexagonal Z(33), but it is of interest that 
Thibault (1943) reports that isotropic carborundum generally contains 
lamelle of «-SiC. It requires something more to explain the absence 
of Z((33), 31), Z((33), 36), Z((33), 35) and of Z((33), 34) apart from 
the one example n=0. The explanation may be connected with growth 
rate as mentioned below, or may be connected with the severity of misfit 
at the slipped surface. It is alternatively possible that the initial 
structure contains a stacking fault in the region whose edge becomes 
exposed : this also gives rise to a polytype with period equal to the step 
height. 

Before the observations of Forty and Verma were available, the 
author sought for an explanation of dislocations predominantly of 
(33), type in the commencement of SiC growth on a quartz crystal (which 
has Si atoms at almost exactly the same separation). The attempt to 
explain the large screw pitches in this way was not successful. 

Whatever the origin of the dislocations, their nature determines the 
crystal structure which can grow on each particular crystal seed. These 
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structures will differ in thermodynamic potential and in growth rate 
under particular conditions of supersaturation. Unstable structures will 
not grow if the supersaturation is small: but probably most of the 
observed polytypes of carborundum differ in thermodynamic potential 
by an amount negligible compared with the supersaturation actually 
occurring. On the other hand, the non-occurrence of structures with 
a 1 in the Zhdanov symbol can well be due to thermodynamic selection. 
The same may apply to those with 5 or 6 in the symbol. 

It remains possible that the cubic structure is actually the thermo- 
dynamically stable one. Since its close-packed faces will form a closed 
body (an octahedrom) it will not normally develop a thin platelet but 
will grow slowly in all directions. Hence, even though stable, it is at 
a disadvantage in growth rate when the medium is also supersaturated 
with respect to the alternative structures. All of these are bounded by 
a parallel pair of close-packed faces only, so that they can grow into 
large thin plates without benefit of dislocations. Their growth rates in 
this direction will be correlated with thermodynamic stability, and the 
evidence suggests that Z(33) is the preferred one of those structures. 


$7. SEPARATED DISLOCATIONS. 


The step generated by slip will usually not terminate abruptly, but 
macroscopically speaking, taper away: on the atomic scale, however, 
the tapering will proceed by discrete dislocations of a whole number 
of SiC layers, and usually a whole period of the initial lattice—six 
SiC layers if this is Z(33). Then the end of the slipped zone consists of 
a number of dislocations, all right-handed or all left-handed and usually 
all of the same kind. Since like dislocations repel each other, they will 
tend to spread apart, out of the line in which they are formed, as growth 
proceeds. This will account for Verma’s observation that in any one 
area the dislocations are commonly all of the same kind. Since this is 
observed both on Z(33) (type II.) and Z(32) (type I.) carborundum crystals, 
we must infer that Z(32), as well as Z(33), is sometimes the initial structure. 
It could, however, be formed secondarily : 7. e. a first buckling and slip 
in Z(33) could generate Z(32) and a second buckling and slip in this 
material generate a number of Z(32) dislocations. 

There is a critical separation of dislocations, of the same order of 
magnitude as the diameter of the critical nucleus of the Gibbs-Volmer- 
Becker—Doring surface nucleation theory, and the spacing of turns of 
a single-branched growth spiral, such that dislocations closer than this 
cooperate in crystal growth, dislocations further apart do not. When they 
cooperate, the growth layers emanating from each of them interleave”. 


*Verma’s fig. 12 illustrates the phenomenon of cooperation. The three 
central dislocations are at such a distance apart that they do not cooperate 
under the high supersaturation of the final stage of growth. In the earlier 
phase, with wider spacing of the spirals, the two above cooperated, dominating 
the one below, as can be seen from the shape of the first widely spaced spiral 
turns. 


1020 F. C. Frank on the Growth 


Three cooperating dislocations, Z(33), Z(33) and Z(32), will generate 
a crystal of structure Z(33 33 32). It will, however, have a misfit surface 
on the line joining the dislocations. The surface energy of this surface 
will oppose their mutual repulsions and, as growth proceeds, gather them 
together. Finally, probably, they will unite in a single hollow core. 

On the other hand, if one of the three dislocations is too far away to 
cooperate, and is dominated by the other two, it will constitute an 
imperfect dislocation of the resulting structure, attached to a fault 
surface extending outwards and tending to draw it away. For example, 
if a Z(33) and a Z(32) dislocation cooperate, the structure will be 
Z(33 32) and the other Z(33) dislocation will bound a fault surface 
extending more or less radially away from the centre of the growth hill. 
If the two Z(33) dislocations cooperate, the structure will be Z(33) and 
the Z(32) dislocation will mark the edge of a fault surface extending 
outwards in the same way. Changes in supersaturation can then lead to 
a change in the structure of the crystal which grows on a given dislocation 
system. But with the progress of growth the dislocations will tend to 
rearrange themselves in a way which reduces the sensitivity to such 
changes. 

The fault surfaces mentioned may act as barriers to the progress of 
growth layers. It is possible that the discontinuity lines seen on Verma’s 
photographs, or some of them, are of this nature. 

The commonly observed coalescence or “‘syntaxis’’ (Baumhauer 
(1915), Ungemach (1935), Thibault (1943), Honjo (1950)) of various 
types of carborundum is now entirely comprehensible. 


§8. CROSS-LACED GROWTH PATTERNS. 

The striking cross-laced growth patterns shown in Verma’s 
figs. 13, 14 and 15 have a simple interpretation. In all these crystals, 
hexagonal or otherwise, the symmetry for any growth layer advancing 
over another is trigonal. For each successive monolayer in the repetition 
sequence of the crystal the force field is slightly different. In a given 
orientation of the step, one monolayer in the repetition sequence will 
have the slowest rate of growth: the rest should therefore overtake it 
and pile up behind it, forming a step of height corresponding to the whole 
repetition period. Now, if the structure is Z(33), 7. e., 
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whichever layer is the slowest growing one in a particular orientation— 
say, the first 4 after a series of Vs—then the third layer after it is 
necessarily the slowest growing one in an orientation inclined at 60° to 
the first. Hence on alternate sides of the hexagon the growth layers 
form alternate groupings—say(VVV 444) and (A44VVV ): 
The growth steps are six layers high (15 A) on the sides of the hexagons 
and three layers high on the zig-zags at the edges of the hexagonal 
pyramids. Fig. 14, showing alternate trigonal spirals, assists the 
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interpretation of figs: 13 and 15, but its explanation is slightly more 
complicated : probably slight surface oxidation has caused the layers to 
retract, separating them into groups only three layers high. 

The conditions necessary for this phenomenon could in principle 
occur on non-hexagonal structures, e.g. Z(32), but are not necessarily 
present as they are for the hexagonal structures. 

Fig. 7, an apparently perfect hexagonal spiral on the hexagonal 
crystal, is actually more difficult to understand. We normally expect 
the rate of advance of a multiple step to be controlled by the deposition 
rate at the bottom of the step: and the symmetry here is trigonal. 
However, the rather high visibility of the phase-contrast image and the 
large central hollow at the dislocation suggests that we probably have 
a dislocation whose strength is a multiple of the repetition period. If 
so, the edge of the step will not be a close-packed surface, but will 
contain re-entrant angles for all orientations. The bottom of the step 
is then not a privileged position and hexagonal symmetry of the growth 
spiral is comprehensible. Indeed, it is slightly under-privileged, with 
respect to diffusion of molecules from the gas. It is this fact which can 
make the step stable against dissociation into steps of height corresponding 
to one repetition period, in a spiral of more than one branch. 
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CIII. The Correlation Function in the Analysis of Directive Wave 
Propagation. 


By Henrik NoptvEDtT, 
Cavendish Laboratory*. 


[Received June 5, 1951.] 


SuMMARY. 

It is shown that the auto-correlation function might be applied when 
investigating the directive properties of waves which cover a wide 
frequency spectrum. Examples are given which clearly show that the 
subsidiary maxima in the calculated polar diagrams are reduced when 
widening the frequency band of the received signals. These subsidiary 
maxima are, further, smaller than those obtained when using ordinary 
phase comparison methods under the same conditions. 

Some possible implications in the field of audio-acoustics are mentioned, 
dealing with the sense of directivity of the human ears. 


§1. [Ivrropuction. 

Durtne the last decade a considerable amount of literature dealing with 
directive wave-propagation has appeared. Ramsay (1946, 1947, 1948), 
Schelkunoff (1943), Silver (1949) and Stenzel (1939) might be taken as 
representative of the methods used and of the results obtained. In far 
the most cases the waves considered are of a single frequency, and the 
direction of propagation is determined by comparing the phase of the 
wave at two or more points in space. The received waves are added 
and, due to the interference, the receiving system will have a definite 
angular discrimination determined by the position of the receiving 
points in addition to whatever amplitude or phase modifications are 
added afterwards. The resulting so-called polar diagrams might be of 
varying shapes and sharpness. Commonly, however, one is interested in 
getting a maximum sensitivity in one particular direction, and in reducing 
the sensitivity in other directions, particularly the so-called secondary 
lobes or subsidiary maxima. 

The single frequency approach is largely due to the importance attached 
to the directive properties of aerials in high frequency and radar work. 
Since all the energy transmitted in these cases is contained in a very narrow 
frequency band, a receiving equipment with the same narrow bandwidth 
is essential in order to achieve optimum signal to noise performance. 

If, however, one is dealing with electro-magnetic fields of the type 
caused by radio stars, or with an audio-acoustic field, the energy in the 
field is distributed over a wide frequency spectrum. If it is required to 
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determine the direction of such fields, it is, of course, possible by suitable 
filtering to suppress all but a narrow frequency band, and then treat 
the problem with single frequency methods. In this case, however, 
the energy contained in the rest of the spectrum is wasted, and the 
sensitivity of the system is consequently reduced. 

It is possible in certain cases to adapt the additive phase comparison 
method to directivity calculations of wideband frequency fields (Stenzel 
1939, p. 27), but it seems natural when dealing with such fields to explore 
the possibilities of the correlation methods. 


§2. PRINCIPLE. 


If a function of time f(t) is multiplied by the value of the function 
7 seconds later, the product integrated with respect to time and averaged 
one obtains the so-called auto-correlation function 


. 
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This function has achieved great importance in communication theory, 
and its characteristics have been studied for various time functions. 
The parameter in the integrand, indicating the delay between the two: 
factors in the product, is in communication analysis altered by suitable 
delay networks, and the measurement of the correlation function is then 
carried out by electronic multiplication in connection with an integrating 

circuit. 

Dealing with functions representing the propagation of waves in a field, 
the parameter can be given a special significance. Considering two 
non-directional receivers in a field (electro-magnetic or acoustic), the 
difference between the times of arrival of the waves at the two receivers 
is given by :— 

MS 2a sin y 
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where 2a= distance between receivers, c= velocity of propagation. 
Realizing the simple connection between the auto-correlation parameter 
and the angle of approach of the incoming waves, it is obvious that the 
correlation function of the received signals at the two receivers is a 
measure of the directivity of the system. By correlating the two signals 
according to equation (1), a multiplicative comparison is substituted 
for’ the additive comparison used in ordinary single frequency phase 
methods. Using the correlation method the phase composition of the 
incoming waves is of no importance. This will be seen more easily if 
equation (1) is transformed according to Wiener’s theorem (Wiener 1930). 
According to this the correlation function can be written as :-— 


bun)= | (GCOS G27, BG)y) rene se aes (2) 


where ¢(w) is the frequency power spectrum of the incoming wave. 
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Since the angle of approach y is contained in the parameter 7, a 
calculation of ¢,(7) will give a polar diagram of the receiving system, 
depending upon the frequency characteristic of the incoming waves and 
upon the frequency response of the receiving system. For the sake of 
simplicity, it will be assumed in the following that the sensitivity of the 
receivers is constant over the frequency band in question and zero outside. 
Likewise, it will be assumed that the energy distribution of the incoming 
waves will be constant over this band. This condition might be 
approached by suitable filtering. 

If, as was assumed earlier, the two receivers are non-directional, 
¢(w) will be constant over the region of integration, and the calculation 
can easily be carried out. However, if receivers with finite linear 
dimensions are introduced, they will have a directivity of their own 
depending upon the frequency of the waves received. The frequency 
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response of the receiving system for a certain angle of approach y will 
then be determined by the directive properties of the individual 
receivers. (w) wilt consequently express these directive properties as 
the frequency is varied, and the variation of 4(w) will, of course, depend 
on the angle of approach y. This angle will, therefore, appear as a 
parameter in ¢(w) as it appears in the other term in the integrand. 

As will be known, calculation of the directive properties of electro- 
magnetic and acoustic radiators or receivers leads to integrals of the 
following type (Silver 1949, chap. IIT. etc., Rayleigh 1896) : : 


N= | n(ay)* 
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In the electromagnetic case N expresses one of the field components 
at a point in space, (xy) the value of the component at the surface. 
In the acoustic case N is the value of the velocity potential, while 
n(wy) expresses the normal velocity amplitude distribution on the 
vibrating surface. 
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Qa 
k= at —=wave constant. 


r=distance from integrating element to point in space. 


F=radiating or receiving surface. Multiplicative constants have 
been ignored. 


Since, however, ¢(w) in equation (2) involves the power spectrum of 
‘the incoming signals, the expression in equation (3) has to be squared 
before inserting. The expression for the correlation function then takes 
the following form: 


by(7)= | Ni vonage ero: 2 oye Se (4) 

It will be assumed that the distance between the receivers and the 
source of the radiation is great compared with the linear dimensions of 
the receivers, and also great compared with the distance between the 
receivers, such that the incoming waves might be considered as plane. 
‘This is closely approximated in most practical cases. The solution of the 
integral in equation (3) has with these limitations been obtained for a 
variety of surfaces, and the solution of equation (4) is consequently 
simplified. Before dealing with directive receivers, however, the non- 
directive case will be considered. 


§3. APPLICATIONS. 


Case 1. Two non-directional receivers. 
In this case the integration of equation (4) gives 
sin Mwr 
Awrt 


where w,)=centre frequency of received bandwidth and 24w=received 
bandwidth. 

It will be seen from equation (5) that the directivity of the system, 
apart from being determined by the distance between the receivers, is 
also influenced by the bandwidth. As was to be expected, the strength 
of the received signal is proportional to the bandwidth, but in addition, 
the shape of the polar diagram is modified. On fig. 2 the value of 
(sin Jw7/ Mw) cos@_7 is plotted against w 97 for various values of dw. 
The first subsidiary maximum in the polar diagram is clearly reduced by 
increasing the bandwidth, and will for the limiting case of dw=w, take 
the value 21-6 per cent. It is also seen that a considerable bandwidth 
is necessary to get a noticeable reduction. The first subsidiary 
maximum in per cent of the main lobe sensitivity is plotted in fig. 3. 
Even for a reception of an octave of the spectrum, the subsidiary 
maximum is only reduced to 84 per cent of the single frequency 
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value. Since the curves for the primary maxima do not vary appreciably 
with varying bandwidth, this part of the function is plotted for the: 
single frequency case only. 
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Directivity functions for two non-direc- Reduction of secondary lobes with in- 
tional receivers. creasing bandwidth. ° 


Comparison of phase and correlation methods, 


As a matter of interest, the directivity curve for the two non-directional 
receivers, calculated by Stenzel (1939, p. 28) is plotted in fie. 4. The 
curve is given by the equation ‘ 


1 1sin wp)r 
pu sake 0 
R= /(5 -- 5 ie cos or) K Se ear sO (6), 
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for dw=}w,. The directivity curve based on correlation is also plotted 
for the same bandwidth for comparison. While the correlation curve 
obviously approaches zero for increasing wor the curve given by equation 
(6) approaches $4/2. This limiting value will be independent of the 
bandwidth as long as two receivers are used. 


Case 2. Two directional receivers. 


As was mentioned, the sensitivity at a certain angle of approach of a 
directive receiver will vary with frequency according to equation (3). 
The function 7(ay) expresses how the sensitivity of the receiving surface 
itself varies with space coordinates independent of frequency. 

In evaluating the integral in equation (4) a particular case will be taken 
in which the receiving surfaces are identical and rectangular, and the 
sensitivity function 7(ay) constant over the surfaces. Considering points 
at great distances from the receivers in the plane containing the surface 
normal and one of the sides of the rectangle, the directivity in this plane 
is given as 


sin x 
R= ; 
. x . 
where 
w ; 
r= — 2b sin y=ar’, 
c 


2b=side of rectangle, 
y=angle between surface normal and direction of approach. 


Squaring this expression and inserting in equation (4) the following 
integral has to be solved: 


We i \ 2 
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The solution of this is straightforward but somewhat tedious. Inte- 
grating by parts one obtains 


A sin wT\2 a (%@7 sin x Fre 2—a eee sin x ie 
= a a ine ae a7 - 
men OTe eee 2wT J 9 0 
2 Oo (2+a)wr sin x Wz ; 
27 Jo a8 , 


Here is 
w5—upper frequency of received bandwidth. 
w,—lower frequency of received bandwidth. 
4=7/7'=ratio of distance between the centre lines of the receivers: 
to the linear dimensions of the individual receiver. 


The integral might also be solved by using the convolution theorem.. 
A rectangular gating function of width (w,—w,) is introduced in the 
integrand, and the limits are increased to infinity. The Fourier transform 


4A2 
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of the function (sin a/z)? is a triangle, while the transform of the gating 
function will be an expression of the form sina/x. Convolving these 
transforms in the time domain leads to the solution. 

Choosing various values for « the polar diagrams (directivity functions) 
for the receiving system are plotted on figs. 5, 6, 7, using the tabulated 
values for the Sine Integral and reducing the maximum amplitudes to 
unity. The ratios between the upper and lower frequencies of the 
band are in all cases given the values 1, 3, 10, ©, corresponding to 
bandwidths of 0, 50, 91 and 100 per cent. 

The reduction of the subsidiary maxima is most predominant when 
the two receivers are close together, 7. e. when « is small. As will be seen, 
the secondary lobes might almost be brought to disappear when «=1 


Fig. 5. Fig. 6. 
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or «=2, by a sufficient widening of the spectrum. For big values of « 
although the relative reduction is quite considerable, the value of the 
secondary lobe is still a high percentage of the primary maximum, as is 
shown on fig. 7 for «=10, In fact, for increasing distance between the 
receivers, the secondary lobe is approaching 21-6 per cent of the primary 
when widening the bandwidth to the maximum possible extent, as in 
the earlier described case. This is of course due to the fact that the 
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fairly slow variation of the directivity of the individual receivers is 
negligible compared with the rapidly varying term caused by the separa- 
tion of the receivers. © 

A plotting of polar diagrams for the receiving system for values of « 
greater than 10 is therefore not considered necessary, unless the other 
subsidiary maxima, where the individual directivity makes itself felt, 


have to be investigated. 


Fig. 7. 
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It ought to be noticed that even in single frequency fields the polar 
diagram based on correlation has smaller secondary lobes than that 
based on phase comparison. In the case of «=1, the first subsidiary 
maximum would in the latter case be 21-6 per cent. From fig. 5 it is 
realised that this is considerably more than the value reached by the 
curve for w,/w,—1. The reason for the reduced secondary lobe by using 
correlation in this case, is due to the multiplicative comparison resulting 
in a squaring of the original directivity curve. 
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The preceding calculation may be carried out for receivers of various 
shapes. For two receivers with circular surfaces, the individual directivity 
curve is given by 
J (kr) 

aes Wy 


where J, (A7)—=Bessel Function of first order, 


r=radius of receiving surface. 


This curve differs from sin a/a in being slightly blunter for the same 
argument, and in having smaller subsidiary maxima. It will, however, 
not give radically differing results when introduced in the calculations. 
The calculations get increasingly difficult when receivers with varying 
sensitivity over the surfaces are introduced. Likewise, a correlation of 
more than two received signals will in most cases be too complex to be 
analysed. 


§4. CONCLUSIONS. 


It has been shown that correlation methods can be applied to directivity 
calculations, especially dealing with waves covering a wide frequency 
spectrum. The secondary lobes are reduced by increasing the bandwidth 
and might in certain cases take negligible values. To get a substantial 
reduction of these lobes, however, bandwidths of 50 per cent or more 
must be used. It is therefore doubtful if the method will be of practical 
importance in high frequency problems, even if a reduction of secondary 
lobes is also obtained in single frequency cases. 

In audio-acoustics, however, bandwidths of several octaves are 
frequently encountered. In fact, single frequency fields seldom occur. 
The sensitivity of the human ear to sound waves stretches from 16 c./sec. 
to 16,000 c./sec., for a person with a normal sense of hearing. Over this 
frequency range a pair of human ears has a surprisingly accurate sense 
of direction (Kock 1950). Experiments with single frequency sound 
sources have lead to the belief that the sense of direction is based on 
phase comparison, since in this case the determination of direction 
coincides with the direction of phase equality. It is known, however, 
that the direction of pulses, transient waves and complex noises generally 
is more easily determined than is that of continuous single frequency 
sound, In the case of those complex wide-band frequency noises, phase 
equality will exist for several frequency components for various angles, 
apart from the ahead direction, and one ought to expect that the sense 
of direction would be impaired. Attempts to explain the sense of direction 
at higher frequencies have taken diffraction phenomena around the 
head into account, and this might be the resason for the experienced 
results. 

On the other hand, if the sensitivity of direction of the human ears 
is based on a sort of correlation mechanism, the phase composition of 
the incoming sounds would be of no importance and the sense of direction 
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‘would automatically be improved by an increasing bandwidth. For 
‘single frequency sounds the correlation would lead to the same results 
as experienced in the experiments carried out. 

One further argument for the correlation point of view is the fact that 
-a sound is detected more easily against a “‘ white noise’ background if 
the sound and noise sources are angularly displaced, 7. e. the masking 
-effect of the noise source is gradually decreasing if it is moved away from 
the ahead direction. It is difficult to see how this can be explained on a 
pure phase basis, while according to the preceding treatment the influence 
of the auto-correlation function of the noise received at the two ears will 
‘be reduced when the noise source is angularly displaced from the sound 
source. 

ACKNOWLEDGMENTS. 


This work was done during the tenure of a fellowship from Norges 
‘Teknisk Naturvidenskapelige Forskningsrid, to which institution 
_.acknowledgment is paid. Thanks are also due to Mr. J. A. Ratcliffe, 
‘Cavendish Laboratory for helpful suggestions and discussions. 


REFERENCES. 


Kock, W. E., 1950, J. Acoust. Soc. Amer., 22, 801-804. 

Ramsay, J. F., 1946, Marconi Review, 9, No.4; 1947, Ibid., 10, Nos. 1, 2, 3, 4; 
1948, Ibid., 11, Nos. 1 and 2. 

RayYLeIcH, 1896, Theory of Sound, Vol. II., p. 107. 

\SCHELKUNOFF, S. A., 1943, Bell. Syst. Tech. J., 22, 80-107. 

Stivers, 1949, Microwave Antenna Theory and Design (M.I.T. Radiation Labor- 
atory Series), Vol. 12. 

Srenzev, H., 1939, Leitfaden zur Berechnung von Schallvorgdngen (Berlin : 
Springer). 

“WreneEr, N., 1930, Acta Mathematica, 55, 117-258. 


fealOpeane| 


CIV. Masses and Modes of Decay of Heavy Mesons*.— 
Part I. «-Particles. - 


By C. O’CEALLAIGH f, 
H. H. Wills Physical Laboratory, University of Bristol f. 


[Received July 15, 1951.] 


[Plates XXIX. & XXX.] 


SUMMARY. 


Observations by the photographic method are described which 
establish the existence of heavy charged mesons, of mass about. 
1200m,, which decay with the emission of a single charged particle. 
In one case, it is shown that the secondary charged particle is a y-meson 
of energy 5-9 MeV. In another, it is a particle of charge |e| and of energy 
about 200 MeV. If the two particles are of the same type with similar 
modes of decay, then the charged secondary particle must be accompanied 
by at least two neutral particles. The question whether the «-particles 
are identical with the 7-mesons which decay into three charged particles, 
is discussed. 


§1. Iyrropuction. 


Four events have recently been found in this laboratory, recorded in 
photographic plates exposed to the cosmic radiation, which represent 
the decay of charged particles of mass ~1200m,. The particles were 
at, or very close to, the end of their range when they transformed into 
a single particle of charge |e| and one or more neutral particles. It will 
be convenient to refer to such unstable heavy particles with this mode of 
decay, as x-mesons. 

The four examples found hitherto have all occurred in plates of type 
Ilford G5, with emulsion thickness 400, exposed on the Jungfraujoch 
high altitude station under a lead plate 30 em. thick. Two of them, which 
we shall refer to as x, and «,, respectively, give particularly favourable 
conditions for measuring the mass and energy of the particles, and it is 
with a study of these two particles that this paper is chiefly concerned. 

The first example of these heavy mesons was found while studying the 
distribution in energy of the electrons emitted in the decay of .-mesons. 
Superficially the event appeared to be of this type, but scattering 
observations showed that the energy of the secondary particle, if an 


*The main features of this and the following paper were given at the 
Copenhagen Conference on Quantum Physics on July 5th, 1951. 

t On leave of absence from University College, Cork. 

t Communicated by Professor C. F, Powell, F.RB.S. 
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electron, was ~240 MeV. This is more than four times the maximum 
energy of the decay electrons from u-mesons, and the characteristics of 
the primary particle were therefore examined. Its track was then found 
to provide unambiguous evidence that the primary particle was several 
times more massive than a u-meson, but less massive than a proton. 


§ 2. CHARACTERISTICS OF THE PRIMARY PARTICLE «yj. 

A photo-micrograph of the event is shown in Pl. XXIX. The primary 
particle decayed near the glass of the plate. It had entered the emulsion 
through the surface, and its path could be traced in the emulsion of 
another plate facing it. The track of the secondary particle also passed out 
of the surface of the emulsion, but near one of the cut-edges of the plate 
so that it was lost to further observation. : 

The point of decay of the primary particle is at P in Pl. XXIX. That its 
direction of motion was towards P is proved by the trend in the values. 
of both multiple scattering and grain density. The appearance of the: 
end of the track is that of a particle which had come to rest before decaying, 
and in what follows, it will be assumed to have done so. 

Since the length of the track exceeded 4000, it was possible to obtain 
a good estimate of the particles’ mass by measuring the multiple: 
scattering as a function of residual range. In determining the mass, 
allowance was made for the variation of energy along the track by a 
method developed in this laboratory by Menon and Rochat (1951). The 
track was long enough to allow four independent estimates of the mass. 
of the primary particle to be made. The measurements were repeated 
by three different observers, and gave a mean value of 1260+290m,. 

An independent estimate of the mass of the primary particle was made- 
by measuring the “ gap-length ” in the track as a function of residual 
range (Hodgson 1950), and for this purpose the emulsion was calibrated 
by making gap-counts on the tracks of four identified »-mesons and 
six protons. These particles all travelled from the direction of the surface. 
and, following a similar path, came to rest in the emulsion at about the- 
same depth as the primary particle «,. This precaution was taken to: 
make the comparison as exact as possible. It was necessary because: 
overfixing had produced some “etching ” of the tracks in the surface 
layers of emulsion down to a depth of about 50 py. 

The measurements of gap-length were confined to regions >50y from 
the surface, and the results are shown in fig. 1. Mean curves are given for: 
the total gap-count (in eye-piece scale divisions) as a function of residual 
range, for ~-mesons (u) and protons (P). The curve of total gap-length 
for the primary particle is denoted by x. This evidence alone proves that: 
the decaying particle is more massive than a 7- or a p-meson. The 
estimate of its mass was derived from the results shown in fig. 1 by a 
procedure due to Perkins (1948). Further, an empirical relation, deduced 
from the experimental results for the »-mesons, allowed the mass of any 
particle of unit charge to be expressed in terms of the total gap-count 
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in its track—as measured between any selected pair of values of the 
residual range. This gave a second estimate of the mass. A number of 
independent values were obtained by both procedures and showed 
satisfactory agreement. The mean value was 1385+200m,, and this is 
in accord with the value 1260-290 m, obtained by the independent method 
based on observations of scattering and residual range. Giving equal 
weight to the two results, we adopt a value 1320+170m, for the mass 
of particle «,. 
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$3. THe SECONDARY PARTICLE. 


The visible length of track of the secondary particle was 2200, and 
its grain-density 18-5--0-9 per 50, a value not significantly different 
from the minimum value observed in the same plate; namely, 
Jmin=17-+0-8. The mean deviation (fé), due to multiple scattering, 
was very carefully determined by several observers using different 
microscopes. The results, which were consistent, were analysed using 
‘several cell-sizes. The mean value, reduced to the standard cell-size of 
1004, was (4) y99=0-:10+0-014°. This allows an estimate of the value 
‘of the quantity pB for the particle. We have pBc=K 19/4) 199, Where 
Koo is the so-called “scattering constant” for 100, cells. We set 
Kyjo9=25 giving pBc=250 MeV. (Gottstein et al. 1951). 
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The above values of (4) 49) and grain-density are such that the 
probability that the mass exceeds 400m, is negligible. The scattering 
measurements therefore show that the momentum of the particle was 
within the interval from 250-300 MeV. /c. 


§4. NATURE OF THE NEUTRAL PARTICLES. 


The conservation laws of momentum and energy demand that one or 
more neutral particles be emitted in the decay of the parent particle. 
If a single neutral particle is produced, an approximate value of its mass 
may be determined if the nature of the charged secondary particle is 
assumed. The results are summarized in Table I. It will be seen that the 
uncertainties in the mass of the primary particle, and of the momentum 
-of the secondary, result in large limits of error for the estimated mass 
-of the neutral particle. 


TaBe I. 
Assumed nature of 
charged particle a bs 
ass of neutral particle 430 +370 5A0 aes 


(Me) 


The limits of uncertainty given in Table I. have been estimated by 
assuming that the experimental measurements of m, and  (¢)j99 are 
each in error by one standard deviation, and then choosing the 
-combinations of these which yield extreme values. Thus, we see that if 
a two-body decay is assumed, the experimental results do not exclude 
any of the following modes of decay : 


K,—>7 OF pw +p (1) 
76 eA Ie PAE) 

Vig COONTA) Les epautaiie « 3fte(3) 

>e +y, (4) 


where v, 7) and V, represent, respectively, the neutrino, the neutral 
7-meson and the lighter neutral V-particle, evidence for the existence 
of which has been put forward by Rochester and Butler (1947) and by 
Armenteros et al. (1951). 


§ 5. PARTICLE kg. 

Photo-micrographs of the second example of the decay of a heavy 
meson, «>, are shown in Pl. XXX. Superficially the event appears to 
represent the successive decay of a 7-particle, 7—1-e, but the -meson has 
the exceptional range of 1098 pw. 
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The range of a u-meson emitted in the decay of a 7-meson is very 
sensitive to any small residual velocity of the parent particle, and further, 
it is well known that the departures from the mean of the values of the 
range commonly observed, can be attributed to straggling (Powell 1948). 
The residual range of the 7-meson at its point of decay in a solid material 
must commonly be less than a fraction of a micron. It must be 
anticipated, however, that decay in flight will sometimes be observed in 
the emulsions of a second plate facing it, and at first sight it appears: 
reasonable to attribute the present observation to such a process. It is,. 
however, very difficult to maintain this view. 

The observed range of the »-meson in the present event is 1098p. 
Since its range is much greater than that resulting from the decay of 
a m-meson at rest it must, if due to such a process, have been thrown 
forward, and not backward, with respect to the line of motion of the: 
parent particle. It will be seen, however, that the »-meson was apparently 
ejected backward. To overcome this difficulty it is necessary to make 
the unlikely assumption that immediately before its instant of decay the 
parent particle was scattered backwards through a large angle, 
subsequently decaying before the change in its direction of motion had 
been made manifest on its track. 

In view of these considerations, the track of the parent particle was. 
studied in detail. The track x, is only 350, long, but the particle entered 
at the surface and its path can be traced in the facing plate, its length 
in this emulsion being ~5800,. Determinations of the mass of the 
particle by methods similar to those described for «,, give a mean value 
of 1125+140m,. Hence, in this case also, we are dealing with the decay 
of a particle several times more massive than a 7-meson. 


§6. NATURE OF THE CHARGED SECONDARY PARTICLE. 


It has been assumed that the secondary particle is a u-meson and the 
evidence for this view is very strong. Its track is sufficiently long to: 
allow estimates of its mass to be made, and the value thus obtained is. 
between 200 and 300m,. If it is a particle of which the existence is 
already established, it must therefore be either a 7- or a p-meson. The 
particle emitted at the end of its range has a specific ionization not 
significantly different from the minimum value. The direct B-decay of 
m-mesons has never been observed, and in an extended study, by Menon, 
Muirhead and Rochat (1951), of more than 3000 7--stars of all types, 
only two cases have been found where relativistic particles were emitted 
(electrons ~5 MeV.). These reasons seem sufficient to reject the 
possibility that the parent meson is a 7-particle. On the other hand, the 
observations are completely consistent with the view that the secondary 
particle is a ~-meson which decays into an electron. The track of this: 
final particle is too short, ~150,, to allow further support for this view 
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to be obtained by a measurement of its energy, but the straightness of 
the track suggests that this was. greater than ~10 MeV. 


§7. NATURE OF THE NEUTRAL PARTICLES. 


If a single neutral particle is assumed to provide the momentum and 
energy balance in the transmutation, its mass can be determined with 
much greater precision than in the case of x, In the present event the 
momentum and mass of the charged secondary particle are accurately 
known, and the errors arise almost entirely from the uncertainties in the 
mass of «,. Further, the momentum of the charged particle is so small 
that most of the residual energy must appear as the rest—-mass of the 
secondary neutral particle. The result obtained for the mass of the 
neutral particle is 900+130m,, and the ““Q” value of reaction is 
~7MeV. Thus decay of x, to » accompanied by either v or zy is excluded, 
but it is possible to identify the neutral particle with the neutral V particle 
of mass 800-900m,. Armenteros et al. (1951) describe a remarkable 
case of the decay in flight of a charged V-particle which was apparently 
associated with the decay in flight of a V°-particle. At first sight it would 
seem reasonable to identify the postulated two-body decay of «, with this 
process. However, calculations show that the Q-value of the cloud- 
chamber event must be several times greater than that of the decay 
of xy. It seems very unlikely, therefore, that both events are examples 
of a similar two-body decay of the same type of particle. 

In addition to the events x, and x,, two other examples of decay of a 
heavy meson have been observed, in which the secondary particle produces 
a track of minimum grain-density, but the geometry of the events is such 
that the mass of the decaying particles cannot be accurately measured. 
Tt is, however, possible to show that they have a mass several times that 
of the z-meson, and it is reasonable to suppose that we are dealing with 
events similar in character to x, and x,. In both cases the grain-density 
of the track of the secondary particle was indistinguishable from the 
minimum. 

Hitherto, no consideration has been given to the possibility that the 
two events «, and x, represent similar modes of decay of particles of the 
same type. Since the estimates of the masses of the two particles do 
not differ significantly, this possibility cannot be excluded. If it is correct, 
both particles must decay with the emission of a »-meson and two or more 
neutral particles. 

If the two particles x, and x, are indeed of the same type, it is at first 
sight surprising that in the first two examples observed to decay at rest, 
there is such a wide disparity in the energy of the charged p-meson. 
This difficulty largely disappears when it is considered that an event 
of the type «x, is much less likely to escape observation than one of 
type x,. It is clearly essential to observe other examples of similar 
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particles which provide tracks suitable for measurement, both of the 
primary and secondary particles, before any final decision on these 
questions can be reached. . 


§ 8. RELATION OF THE x- AND 7-PARTICLES. 


The determinations of the mass of the «-particles and the results of 
similar measurements on the r-particles are summarized in Table IV. of 
the following paper. The results suggest that the «-particles are about 
1-3 times more massive than the 7-particles, but the inaccuracies in the 
experimental values are such that the possibility cannot be excluded that 
they represent alternative modes of decay of particles of the same type. 
Further, they may also be of the same type as the unstable charged 
particles of Rochester and Butler which are observed to decay in flight. 
It has been pointed out by Bethe (1951) that this is the most economical 
assumption, in the sense that it calls for the smallest number of different 
types of heavy mesons. We believe, however, that the experimental 
evidence, while not decisive, is against such a view. The most cogent 
reasons are the following : 

Assuming that «,=«x,=«; and that «=7, it follows that the spins of 
7 and « are equal. Further, since 7 is believed to decay into three 
z-particles, it follows that its spin is integral. 

Suppose that « decays with the emission of a u-meson of half-integral 
spin and two neutral particles. Since « is assumed to be of integral spin 
the two neutral particles cannot be similar. One must be half integral 
and the other integral. The simplest mode of decay of «, involving known 
particles only, is then represented :— 

K>p+v-+7, (i.) where v and 7» represent a neutrino and a neutral 
a-particle respectively. 

The most accurate determination of the mass of a heavy meson appears 
to be that of the 7-particle, viz., 966m,, and assuming the identity of 
z and x, this is also the mass of x. In the event «,, however, the scattering 
<)199 Of the secondary particle corresponds to a value of the quantity 
pB equal to 250 MeV./c; if this secondary particle is a .-meson, its energy 
is therefore 182 MeV. Applying the conservation principles to the mode 
of decay (1), and assuming the neutrino to have negligible energy, the mass 
m, is found to be 1160m,. The uncertainties in the value of CO yion 
for the secondary particle are such that we cannot exclude this mode of 
decay for a particle of mass 966m,. Since, however, we have assumed 
little energy to be taken by the neutrino, it appears improbable that the 
mass of the parent meson «, is as low as 966m. 

These considerations do not exclude the possibility that the various. 
charged V and « particles, which emit a single charged secondary particle, 
represent a variety of modes of decay of particles of the same type. 
But it would be surprising to find several different modes, occurring with 
approximately equal probability, such as the limited experimental 
material suggests. 
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heavv meson, «,. reaches the end of its range and decays with the emission of a single 
Vy Deal) g 


particle of charge lel. The track of «, was longer than 4000 microns and measure- 
ments of scattering and grain-density give a value for the mass of the particle 
1320--170m,. The energy of the secondary particle, if a z-meson, is 180 MeV.; if an 


electron, 250 MeV. 
Observer ; Miss A, Stradling ; Photograph by : A. R. Gattiker and Miss P. Dyer. 
To face page 1038 
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A heavy meson, ky, reaches the end of its range at a point P and decays with the 
emission of aj-meson of range 10984. The p-meson stops at the point Q 
and emits an electron. The track of the j.-meson is given in two parts. 


The mass of the particle as given by measurements of scattering and grain- 
density is 1125+140m,. 
” os ¢ 


Observer Mrs. I. Powell. Photograph by A. R. Gattiker and 
Miss P. Dyer. 
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SuMMARY. 

Further evidence which establishes the existence of heavy charged 
mesons which decay into three charged particles is presented. It is 
shown that the three secondary particles are probably z-particles, and 
that the mass of the parent mesons, called 7-particles, is 966+8m.,. 
The frequency of occurrence of «- and r-particles at the Jungfraujoch, 
in comparison with the number of z-particles, is discussed, and it is shown 
that in the nuclear collisions involving protons of energy greater than 
10 BeV., a large fraction of the energy lost in the creation of mesons, 
must appear as x- and 7-particles. The evidence suggests that the life-time 
of the heavy mesons is about 10-* sec. The most favourable conditions 
for the detection of the particles by the photographic method are 
discussed. 


§1. IyTRODUCTION. 

More than two years ago a description was given (Brown et al. 1949) 
of an “event ’”’, recorded in a photographic emulsion, which appeared 
to represent the spontaneous decay of a charged particle with mass 
~1000m, into three mesons, probably z-particles. Photo-micrographs 
‘of the tracks are reproduced in Pl. XX XI. Although the evidence for the 
interpretation was very strong, it was necessary to observe similar events 
before the existence of heavy mesons with this particular mode of decay 
could be taken as established. 

Support was long in coming, but about a year ago, Harding (1950) 
published a description of two observations similar to the first. In each, 
three particles appeared to diverge from a point at the end of the range 
of another particle, and the directions of motion of the three emergent 
particles were co-planar. Further, one of the two new events allowed 
the mass of the parent particle to be estimated, and the value obtained was 
again of the order of 1000m.,. 

Another example with similar characteristics, which gives exceptionally 
favourable conditions for measuring the mass of the primary particle, 
has recently been found in this laboratory. The mass of the particle 
is 969+10m,, assuming that of the 7-meson to be 274m, (Alvarez 1950). 
‘The observations on the two previous events have been reinterpreted, 

*A letter to the Editor on the same subject by P. E..Hodgson is also 


published in this issue, see p. 1060.—Hditor.] 
+ Communicated by Professor C. F. Powell, F.R.S. 
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and they are shown to be consistent with this value. We therefore 
believe that the existence of these particles with their distinctive mode of 
decay can now be taken as established. It will be convenient to refer 
to them as t-mesons. 

In this paper, the main features of the new observations are first 
outlined and compared with the earlier results; the original evidence on 
the nature of the secondary particles is supported and an improved value 
for the mass of the z-particles deduced. The frequency of occurrence 
of the heavy mesons and the best experimental conditions for their 
detection by the photographic method are then discussed. 


§2. NEw OBSERVATIONS. 


A mosaic of photo-micrographs of the new event which was found in 
an Ilford G5 emulsion, 400 thick, exposed for five weeks at the 
Jungfraujoch under 30 cm. of lead, is shown in Pl. XXXII. The direction 
of motion of the particle which produced the track 7 is unambiguous, and the 
particle was at, or very close to, the end of its range at the point P. The 
total length of the track in the emulsion is 20704. The measurements of 
the scattering of the particle as a function of the residual range gave a 
value of 1015+280m, for its mass. Using, in addition, the information 
given by gap-counting, the best value obtained was 1000+180m.,. 

From the point P, three particles with charge |e] diverged, and they all 
left the emulsion before coming to rest. One particle, however, track a, 
moved nearly parallel to the plane of the emulsion, and it entered the 
glass 6-4 mm. from P. It was then moving relatively slowly, as shown 
by the scattering of the track, and its residual range at the point of exit 
from the emulsion was only about 400. Tracks (b) and (c) dip much more 
steeply and, in the emulsion, they are only 1204 and 490 long, 
respectively. These values represent the true lengths after correction 
for dip. 

The tracks, a, b and c are sufficiently long, and sufficiently straight, to 
allow their orientations to be measured accurately. The observations. 
show that the initial directions of motion of the particles were co-planar 
to within the limits of error of the measurements ; in the emulsion, before 
processing, the direction of any one particle was inclined to the plane 
containing those of the other two at an angle of 2°+2°. This feature— 
also observed in the three events previously described—makes it 
reasonable to assume that the charged particles are the only products of 
the decay of the parent meson. 

The great length of the track (a) makes it possible to determine the mass 
of the particle—by observations on the grain-density and the scattering 
parameter, %—over a succession of segments of the track, and the mean 
of the values thus obtained was 285-+20m,. It is therefore almost 
certain that the particle was not a .-meson, and the result strongly suggests 
that it was either a 7-meson or an unknown particle of nearly equal mass- 
We shall assume the former. 
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Results obtained for the mean angle of scattering, x, and grain-density 
for different intervals along the track are given in Table I. The measure- 
ments show conclusively that the particle had lost almost all its initial 
energy before it entered the glass. By comparing the results of Table I. 
with similar data obtained on tracks of y-mesons and protons which 
stopped in the emulsion, we have deduced a value for the residual range, 
r, of the particle as it left the emulsion of 400-2004. The best value for 
the initial energy of the particle is then determined from its range in the 
emulsion, which is R=6-8+0-:2 mm. The value of the initial energy thus 
determined is insensitive to the errors in the estimate of r; the uncertainty 
in the energy of the particle at P, introduced by these errors, is of the 
same order of magnitude as that due to straggling. Using the range- 
energy relation for protons (Lattes et al. 1947, Bradner et al. 1950), the 
energy of a particle of mass 274m, and range 6-80 mm. can be shown 
to be 19:-0-0-4 MeV.; and its momentum 75-4 MeV./c. 


TABL Ml 


Variation of the grain-density and the scattering parameter, x, along the 
track of particle a, see Pl. XXIX. The ranges referred to are the 
distances measured along the track from the point at which the 
particle entered the glass. 


Interval in 
range 


grain-density 
/50. 
6375u—2575y 4500 0-85 +0-076 51 
261 5u—1335. 2000. 1-26 40-189 63 
500p— 0 250 3-0 +0-45 black 


R a, 
Mean range & 190u 


The relative directions of motion of the three secondary particles in 
their common plane at the time of their formation can be easily determined. 
Knowing the momentum of particle a, and assuming the vector sum of the 
momenta of the three particles to be zero, the momenta of particles 
(b) and (c) can be calculated. The values thus obtained are 
85-8-+1 MeV./c and 98-3+1 MeV./c, respectively. These observations, 
together with the grain-density in the tracks, allow estimates of the masses 
of the particles to be made. It is thus found that the mass of particle (b) is 
240+30m,, and of particle (c), 280+15m,, the limits of error 
corresponding to the standard deviations in the number of grains counted. 
We believe that for particle (b), the observed mass-value is too low owing 
to a tendency to underestimate the grain-density in steeply dipping tracks. 

These observations give strong support for the original assumption 
that the secondary particles are all z-mesons. If this is so, the energy of 
the particles can be deduced from their momenta and the values thus 
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‘obtained are :—particle (a), 19:0+0-4 MeV.; (b), 24-2-+2-0 MeV.: and 
(c), 32-0+2-0 MeV. The total kinetic energy of the three perieies is 
‘therefore 75-2-++5-0 MeV. 


§3. THe Enercy ReLeasé IN THE Decay OF THE 7-PARTICLE. 

In order to compare the present determination of the release of kinetic 
energy with that in the two previous events on which measurements were 
possible, the old observations have been re-examined using the latest 
values for the mass of the 7-meson, m,—274--2m, (Alvarez 1950). 
‘The results are summarized in Table IT. 


TaBLeE II. 
Total Kinetic Energy of the Decay Particles. 


Assumed to be z-particles 


Event 
Particle (a) Particle (b) Particle (c) Total 


] 1:04+0-10 MeV. | 31+4 MeV. 33+4 MeV. 65+8 MeV. 
2 50+7:5 MeV. 13+2 MeV. 22+3 MeV. 
3" 


85 +15 MeV. 
19-+0-4MeV. |24-242MsV. | 3242MeV. | 75-5 MeV. 


c2 
os 
+ 
HS 
= 
= 
a>) 
Be 
$ 


Weighted Mean 7: 


The corresponding value of m, is 966-48 m,., assuming m,—274 m,. 


The method employed for the analysis of event No. 1 was similar to 
that used above, and depended upon the knowledge of the energy of 
particle (a). In event (1) this particle—almost certainly a 7~-particle—is 
of short range, and the estimated value of its energy is subject to errors 
of +6 per cent owing to straggling. Further, tracks 6 and ¢ are inclined 
‘to one another at an angle of only ~10°, and the values of the momenta 
of the corresponding particles, deduced from that of (a) by assuming 

a conservation of momentum, are sensitive to errors in the precise value 
.of this angle. 

‘In event number (2), the determination of the energy depended on 
measuring the scattering parameter for one of the tracks, of length 
1300, and the values thus obtained were subject to greater statistical 
fluctuations than those based on range measurements. 

The weighted mean of the three determinations is 73:-5--4 MeV. 
Assuming all the product particles are 7-mesons the mass of the parent 
particle, m,, is 966-+-8m,. Events (2) and (3) were isolated for detailed 
study because of their resemblance to event number (1); and only after 
a detailed analysis was it established that, in each case, the directions 
of motion of the three secondary particles were co-planar, and that the 
magnitude of the energy release was closely similar to that in the first 
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event. No other event has been found in this laboratory with even 
superficially similar characteristics. It can hardly be doubted, therefore, 
that the events represent similar modes of decay of particles of the same 
type. Otherwise very improbable and arbitrary assumptions have to be 
introduced in order to account for the events in any other way. 


§4. Toe NATURE OF THE SECONDARY PARTICLES. 

The above value of m. depends on the correctness of the assumptions, 
(1), that the three charged particles produced by the decay of the 
r-particle are not accompanied by any form of neutral radiation ; and, 
(2), that they are all z-mesons. The evidence for the nature of the 
secondary particles derived from a study of all three events can be 
presented in the following alternative manner :— 

We have seen that, for any one event, it is generally possible to measure: 
the angles between the directions of motion of the three secondary 
particles in their common plane, and the grain-density of their tracks. 
This makes it possible to determine the relative values of the masses for 
the three secondary particles in each event. Thus, if the assumption of 
a momentum balance is correct, the momentum of one particle is 
proportional to the sine of the angle between the directions of motion 
of the other two. In addition, the grain density in a track gives a measure: 

3 momentum . 
of the quantity B/4/(1—f?)= ees x f(g). It follows that the ratio of 
the rest masses can be written : 


; oes sin oe : sin Pac : sin Pav 

RCA MAY JACK 
where ¢,,, represents the angle between the directions of motion of 
particles b and ¢, etc., and g, is the grain density of track a, ete. 

In event No. 1, only two tracks can be compared in this way, but the 
other, track (a) is due to a 7~-particle of very short range. For tracks 
b and c it is found that m, : m,. :: 1 : 1:02+0-10. 

In event No. 2 (Harding 1950), using data given by the author, we- 
find : 


Mg tM, 1.m,. 2 > Tob 10a TO: 


a iC 
and for the present event, the corresponding results are : 
Ma: M,:m,:: 1: 0-88+0-10 : 1-:03-+0-05. 

The results are shown in Table III. which summarizes information: 
about the masses of all the nine secondary particles on which measurements 
have been made. 

From Table IIT. it will be seen that all the particles on which measure- 
ments have been made give results which are consistent with a unique 
mass, and in three cases the evidence is in fact very strong that these- 
particles are a-mesons. Whilst the possibility cannot be excluded 
that one of the three product particles produced in this mode of decay 
of a 7-particle is a ~-meson, the evidence strongly favours the assumption 
that they are all 7-particles. 
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The co-planarity of the three secondary particles appears to be an 
invariable feature of the transmutation. We cannot, however, exclude 
the possibility that neutrinos, or y-radiation, of relatively low energy and 
momentum, are emitted in the transmutation. Thus, a neutrino of 


TABLE III. 


Absolute or Relative Values of Mass of Secondary Particles. 


Event Particle (a) Particle (b) Particle (c) 
1 7=274m, 280 +30 1-02+12% 
2 1 1:104+15% 1-:0+10% 
3 285-+420m, 0-88+11% | 1-034 5% 


energy ~10 MeV. would produce an average departure from co-planarity 
of the tracks of the three charged particles of only ~3°, and the precision 
of the angular measurements does not yet allow us to exclude it as a 
possibility. If such unobserved particles are emitted, the mass m_ might 
be too low by about 20m,. It appears very improbable, however, that 
in the three observed events which correspond to the decay of the 
7-meson, the assumed fourth particle of a four-body decay should in 
every case have been of exceptionally low energy and momentum. 

The results obtained for the mass of the 7- and «-mesons, using the 
various methods, are summarized in Table IV. 


TABLE LY. 


Mass of Primary Particle in 77g 


Conser- 

Tet vation Mode 
| Particle | Authors rae - avR g.d.uR of of 
or trae momen-| decay 

tum 


T Brown et al} 3100p 990 +270 | 1080 +160 | 949 +16 | 7>7-+-7-+4-7 
T Harding 740 — a= 989 +30 
T Fowler et al| 20704 |1015+280] 9104220} 969+10 


Ky O’Ceallaigh | 4100u | 1260 +290 | 1350 +180 
iy O’Ceallaigh | 56704 | 1125+260] 1125-+-150 Ket 


§5. FREQUENCY OF OCCURRENCE OF T- AND «-PARTICLES. 


Although the statistical weight of the observations is at present very 
low, they suggest that in energetic nuclear collisions a considerable part 
of the available energy appears in the form of heavy mesons. In plates 
exposed under 30 cm, of lead on the Jungfraujoch, five 7- or «-particles 
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(O’Ceallaigh 1951) have been found in a volume of emulsion in which 
750 z-particles are recorded. We may therefore write N,, ,/N,= 1/150, 
where N,. is the sum of the number of 7- and «-particles. For the: 
following reasons, the order of magnitude of this ratio is very suggestive. 

Recent experiments in this laboratory (Camerini et al. 1951) show 
that in the disintegrations produced by protons of energy, E, in the range 
from 2 to 10 BeV., the total energy appearing in the form of emitted 
7-mesons is approximately proportional to E. Let us assume, first, that 
in the nuclear collisions involving protons of energy above 10 BeV., the 
same fraction of the energy appears as mesons; secondly, that this energy 
is equally divided among 7-particles on the one hand, and 7r- or «-particles 
on the other; and thirdly, that the 7- or «-mesons are produced in the 
centre-of-mass system of the interacting nucleons with a distribution in 
velocity and angle similar to that determined for the 7-particles. With 
these assumptions, and knowing the distribution in energy of the fast 
protons and neutrons in the cosmic radiation at the Jungfraujoch, the: 
expected ratio of the numbers of 7-particles produced to that of heavy 
mesons can be calculated. The result is N,./N,~1/75, which is of the: 
same order of magnitude as that observed. 

The above result is particularly striking because no account has been 
taken of the difficulties of recognizing the tracks of «-mesons. These 
particles have to be distinguished from p-mesons which are about 
a thousand times more numerous, and they can only be identified with 
confidence if the track of the primary particle is longer than 2 or 3 mm..,. 
so that the mass can be determined ; or if the secondary particle is of 
energy greater than 100 MeV. and produces a track longer than 3 mm. 
These requirements are so stringent that we can be certain that we fail to: 
recognize at least twice as many «-particles in our emulsions as we 
identify. On the other hand, like the z-mesons, a 7-particle should be: 
recognized with a high degree of efficiency. 

The above considerations take no account of the possibility that the 
«- and 7-particles may be the products of decay of heavier mesons; or 
that the «- and r-mesons are short-lived so that they have an appreciable: 
chance of decaying in flight even when moving in solid substances. Both 
these features would tend to reduce the expected value of the ratio: 
N,,,/N,, and this would make the number actually observed even more: 
surprising. The observations therefore give strong support for the view 
that in the nuclear transmutations produced by nucleons of energy above 
10 BeV., a considerable fraction of the energy available in the collision 
is transformed into «- or 7-particles. 

In the preceding paragraphs we have made no distinction between: 
«- and z-particles because we have no precise knowledge of their relative 
frequencies of occurrence. There are, however, reasons for believing 
that the «-particles are more numerous than the 7’s. 

In the Plates, in which the present events are recorded, four x- and one 
z-particle have been found. We have seen that the mode of decay of the 
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7-meson makes it easily recognized. On the other hand, the «-particles. 
can be distinguished from y-mesons only with difficulty and in favourable 
cases. It is therefore possible that there is a large preponderance of 
«-particles, and if this is so the remarks in the following paragraphs on 
the life-time of the particles and the best conditions for their detection, 
would largely apply to the «-particles. 


§6. Lire-Time oF THE HEAvy Mesons. 


The total time taken by a r- or x-particle to traverse its observed track 
in the emulsion before coming to rest can easily be calculated, and these- 
time intervals, for all the x- and 7-particles identified in this laboratory, 
are of the order of 5x10 sec. We cannot, however, draw any 
conclusions from this observation because we do not know how many 
similar particles decayed in flight. 

The considerations of the preceding section give strong support, however, 
for the view that the heavy mesons have a mean life-time not less than 
10~® sec. Thus, assume a life-time of 5x 10-1 sec., and a distribution 
in velocity of the 7- or «-particles postulated in the previous section. 
It is then found that, if moving in copper, only about 5 per cent of the 
particles would be brought to rest before decaying. For a life-time of 
10~® sec., the corresponding result is about 25 per cent. The actual: 
number expected might be considerably less than this if the particles. 
can also be removed through collisions with nuclei. Since we have already 
accounted for the observed number of «,7-particles without assuming 
any loss by these processes, it is reasonable to suppose that the life-time: 
is not less than 10~° sec. 


§7. SIGN OF THE CHARGE OF THE x- AND 7-PARTICLES. 


There is no direct evidence, provided by the present experiments, for 
the sign of the charge of the x- and 7-particles. We may assume 
tentatively, however, that they are produced in nucleon—nucleon collisions. 
and have a strong interaction with nuclei. It is then reasonable to consider 
the possibility that negative heavy mesons, if absorbed by nuclei, will 
produce o-stars similar to those resulting from the nuclear capture of 
m-particles. Because of the greater mass of heavy mesons, such 
disintegrations might be correspondingly more energetic, and lead to the 
emission of more charged particles than are commonly found in the 
o-stars. No such examples which could be attributed to the nuclear 
capture of a heavy meson have been found in this laboratory, but for the 
following reasons, this result is not decisive :— 7 

There is, first, the possibility that following capture by nuclei, most of 
the rest-energy of the heavy mesons may escape in the form of neutral 
particles or y-radiation. Such is the case in the nuclear disintegrations 
produced by ,--particles. If this is so, the o-stars produced by heavy 
mesons might easily be confused with those due to the capture of negative 
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x-particles. A search was therefore made for such events, in plates exposed 
to the cosmic radiation, by measuring the masses of the particles producing 
the o-stars. In 100 examples examined hitherto none has been found in 
which the parent particle could be identified as a heavy meson. 

It must also be borne in mind that since the x- and 7-particles are much 
more massive than z-mesons the scattering of the particles and the 
variation of grain-density with range are much less than for the 7-mesons. 
This makes it more difficult for an observer to recognize the parent track 
of such an event. Particularly if it is short, he may mistake it for the 
track of an outgoing particle, and regard the event as an ordinary 
cosmic-ray star. 

In addition to these considerations, it is not certain that negative 
7- and «-particles, even if directly produced in nucleon-nucleon collisions, 
will, when captured by an atom, lead to the disintegration of its nucleus. 
The observed events may therefore represent the spontaneous decay of 
both positive and negative heavy mesons. Since the total energy of the 
three 7-particles is not extremely great, and since their directions of 
emission are co-planar, it may be possible, in favourable cases and in 
thick emulsions, to identify all the product 7-mesons, and thus the sign 
of the charge of the parent particle. 


- $8. DresigN oF EXPERIMENTS FOR DETECTION OF x- and 7-PARTICLES. 

The frequency of occurrence of heavy mesons in nuclear disintegrations 
of great energy, and their significance for the development of meson 
theory, make it a matter of great importance to observe many of the 
particles so that their properties may be established. They occur so rarely 
that such a programme would seem to demand, for its successful 
completion, the joint efforts of experimenters all over the world using all 
available technical devices. It is therefore appropriate to consider the 
most suitable design of experiments for the recording of the particles 
by the photographic method. 

The principal difficulty which has to be met in studying the «- and 
t-particles is that they must be detected amongst a much greater 
number of particles of different types. Thus in searching for «-particles 
it is desirable to examine the tracks of all protons and p-mesons which 
stop in the emulsion. This work would be more fruitful if the proportion 
of «-particles could be increased. Since the heavy mesons can only be 
produced by particles of great energy, and since disintegrations by 
particles of energy below ~5 BeV. only contribute to an undesirable 
“ background ” of tracks, the most favourable experimental conditions 
would be provided by experiments at great altitudes at geomagnetic 
latitudes less than 40°. At such latitudes, the ‘‘ cut-off ’’ ensures that the 
minimum energy of the protons reaching the top of the atmosphere is 
about 5 BeV.; and by going to great altitudes we ensure that nuclear 
collisions in the atmosphere have not produced a large flux of particles 
of lower energy. 
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heavy meson, 7, comes to the end of its range in the emulsion and decays into 3 charged 
particles, one of which is az~-particle of short range which produces a nuclear disinte- 
gration. There is strong evidence that the tracks b and ¢ were also produced by 
q7-mesons. The direction of motion of the three secondary particles are co-planar. 


Observer Mrs. Van der Merwe. Photograph by A. R. Gattiker and Miss P. Dyer. 
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A heavy meson, 7, reaches the end of its range at the point P and decays with tl 
emission of three charged particles, probably z-particles. The grain-density a1 
scattering in the track of the particle, 7, gives a measure of its mass and the val 
thus found is about 1000m,. The tracks of the secondary particles are co-planar. 


Observer Mrs, IT. Powell. Photograph by A. R. Gattiker and Miss P. Dye 
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In such experiments on heavy mesons it is necessary to surround the 
photographic plates by local matter in which the particles can be generated. 
The best design for this assembly depends on the lifetime of the particles 
-and this, it has been suggested, is greater than 510-1 sec. Further, 
the experiments with expansion chambers on the charged V-particles, 
in which they are observed to decay in flight, suggests that the lifetime 
is not much greater than 10-® sec. Taking this as the best value at present 
available, it may be shown that the most favourable conditions for 
observing the heavy mesons are given by surrounding the photographic 
plates by material of medium atomic weight, such as copper, to a thickness 
-of about 10 cm. Such an assembly would weigh about 80 lb., and it 
could be carried to great altitudes by balloons. There would be some 
advantage in increasing the thickness of the surrounding material, but 
it would be offset by the great increase in the load, and the reduction 
‘in the altitudes which could then be reached. 


§10. CONCLUSIONS. 


The most important conclusions which follow from the present 
‘experiments may be summarized in the following terms :— 

(a) Taken in conjunction with those of Harding, the observations 
confirm the existence of charged mesons of mass ~1000m,, called 
t-particles, which decay into three particles of charge |e. 

(b) Most of the secondary particles are 7-mesons and the evidence is 
very strong, but not finally decisive, that they are all of this type. 

(c) If the 7-meson decays into three z-particles its mass is 966-+-8m,. 

(d) The lifetime of the heavy mesons («,7) is of the order of 10~° sec., 
and good conditions for their detection, using the photographie 
method, can be obtained by making high altitude balloon flights 
at geomagnetic latitudes less than 40°. 
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Distributions of Second Differences and the Technique of Overlapping 
Cells. 
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SUMMARY. 


Predictions based on a theoretical treatment (Moliere 1951) of the 
problem of overlapping cells, are in satisfactory agreement with the 
results of experiment. They would indicate that the dispersions of the 
sampling distributions of certain statistics commonly used to estimate 
the scattering of particles, are almost independent of A, the coefficient 
of overlap. The dispersions are found to differ but little from those of 
measurements based on a smaller number of non-overlapping cells which 
are shown to be independent in the sense of ordinary statistical theory 
of small samples. It is concluded that the advantage of employing the 
technique of overlap is inappreciable. 


§1. DEFINITION OF OVERLAP. 


To obtain the best value for D,,, the mean absolute value of second 
difference for cell-size ¢ microns, it has been common practice to make use 
of measurements on overlapping cells. This is done by measuring the 
coordinates of the track at intervals c’, a submultiple of c, such that 
c’=c/A. We refer to c’ as the prime cell-size. The overlapping values 
of second difference are then obtained by using all those pairs of prime 
cell coordinate readings which are separated by an interval c. 

Overlap is defined to exist for A>2, so we term A the coefficient of 
overlap, and use it as a convenient measure of the degree of overlap. 
Suppose that %, ¥1, Yo... Y;--. Yy, ave the track coordinates at 
prime-cell interval «’. The following will be the expressions for the 
jth second difference. 


Without overlap, Djx=y,;-;—2y;+Yy;,, (cell-size=c’). si eae, pase ae 
With overlap, Dy, =y;-,;—2y;+y;4,  (cell-size=Ac’). ae Pa cnmoueel ae 


Consider now, a track of length=L and any chosen cell-size c,.. Sequences 
of second difference corresponding to various degrees of overlap may be 
obtained by suitably reducing c,/A, the interval between readings. The 


* On leave of absence from University College, Cork. 
+ Communicated by Professor C. F. Powell, F.R.S. 
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sae of readings of second difference thus obtained will be: 


N= oa —(2A—1). These, however, will not be statistically independent. 


_ Since the labour of computing the differences is increased by a factor 
of order A, it is important to enquire whether the use of overlap will 
result in any commensurate decrease in the dispersion of the experimental. 
results. 


§2. SravTistics oF SaMPLING DisTRIBUTIONS. 


Consider a universe of monoergic particles, and suppose that we may 
select from it tracks all of length=L measured in cell-size=c,. Each 
track regarded as a sample from this universe, will contain n values of 
second difference D;. Let conditions be such that these n values are 
statistically independent. 

As shown in (I.), any measure of the dispersion of the above universe: 
may also be used to measure the energy of the tracks. We choose 
ahy parameter “ p”’ characteristic of the universe which is a convenient 
measure of the energy. The estimate of p derived from the individual. 
samples is called a statistic, and will be distributed about this universal. 
value. The sampling distribution of the chosen statistic is known in 
certain cases. To determine satisfactorily the energy of any individual. 
track, the standard deviation o, of the sampling distribution must be 
small. For the present purposes we consider only the following measures 
of dispersion. For these, except for very small samples, theory shows that 
o, will depend sensibly on n~* where n is the number of independent 
readings in each sample. 

(1) The mean deviation of D,;, from assumed universal mean=0. 
This is the statistic commonly used, and has already been termed D. It 


is defined as 


—— Ip we 
Mine este |), | spe es eens eh 4.) a, (8) 
Nim 


(2) The mean deviation m of D; from the sample mean 
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(3) The standard deviation 
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a 
We distinguish the universal parameters by Greek letters, setting the 
standard deviation=c and the mean deviation=p. The distributions 
of the statistics (4) and (5) are known. That of (1) however, although of | 
greatest interest, appears not to have been studied. In cases where 
slight C-shaped distortion of the emulsion is suspected, it seems best to 
use (4) or (5). In essence this has already been suggested by: 
Fowler (1950). 
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$3. STANDARD DEVIATION OF SAMPLING DISTRIBUTION FOR 
OVERLAPPING CELLS. 

Moliére (1951) has examined the problem of determining the standard 
_deviation of the sampling distribution, when the individual values, of 
which each sample is composed, are not statistically independent. 
Adapting his method to the present case, we find that the second 
differences obtained with overlap are strongly correlated. Table I. gives 
as a function of increasing A, the values of the elements a,;, in the jth 
row of the coupling matrix. 


TABLE I. 


Correlation Coefficients @;,.. 


r l 2 3 4 
k=j—7 ) ) 0 1/256 
j—6 i) 0 0 8/256 
4—5 ) 0 1/108 27/256 
j-4 0 0 8/108 1/4 
j-3 0 1/32 1/4 121/256 
j—2 0) 1/4 60/108 174/256 
q—l 1/4 23/32 93/108 235/156 
j 1 1 l 1 
j+l 1/4 23/32 93/108 235/256 
j+2 i) 1/4 60/108 174/256 
j+3 0 1/3 


/32 1/4 121/256 


For the degree of overlap 4, the general term a;,, k=j—v will be 


2A—v)®—4(A—v)8 
pes fedsey) eet v) = v) 4 v<a 


and 
=A(2A—v)3/403, v>a. _., ia are 


In each sample there will be N, values of second difference. We seek 
to determine the number of these which may be regarded as independent 
from the viewpoint of the theory of sampling. 

Using the methods of large-sample statistics, and assuming that the 
universe is normal, we may deduce from (6) the following expression 
for oy,/o; oy, is the 8.D. of the sampling distribution of s(A); s(A) is 
the statistic “ standard deviation of D, defined by (2) ” 


Oya [| (A) ]1/? 
“ =| san che, Siege CRN Als 


%(A) is a numerical factor which is found to have the following values. 
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N,/A=L/C,—2+1/A, will not differ appreciably from L/C,—1, the 
number of second differences obtained without overlap, and the factor 
(A) decreases very slowly with A. It follows that sharpening of the 
sampling distributions obtained by increasing A is quite negligible. We 
conclude therefore that the advantage of using overlapping cells, is 
not commensurate with the increased labour involved. 

The ratio o,/o for samples of N independent readings from a normal 
universe is known to be (Kenney, (a) Mathematics of Statistics IT.). 


! 2 3 e ON)L/2 : 
6,/o== IN 8h? len? °° a (DAN) a. da (8) 


Comparing (7) with (8) we infer that for an overlap A a fraction 
g(A)=1/As(A) of the second differences may be regarded as independent 
from the viewpoint of the statistical theory of sampling, as the quantity 
N ,g(A) plays the same role in formula (7) as the quantity N in formula (8). 
For the case A=1 this fraction is 1/1-125 or 89 per cent. For values of 
A>1 the rate of variation of g(A) and N, with A is such that their product 
remains practically constant. Thus the standard deviation of the 
sampling distribution will be sensibly independent of A. 


TABLE II. 
Xr ] 2 3 4 
wei teet = A a a a EE  E ————————————————————— ee 
(A) 1-125 1-10 1-076 1-045 


These conclusions depend on two assumptions which in fact are not 
rigorously satisfied. (a) that the universe is normal, and (b) that the 
methods of large sample statistics may be applied to samples of the size. 
commonly used, viz. 10-30 independent cells per sample. Nevertheless, 
we conclude that they are not seriously in error, as will appear from the 
results of the experimental investigation to be described. 


§ 4, EXPERIMENTAL STuDY OF SAMPLING DISTRIBUTIONS 
WITH AND WITHOUT OVERLAP. 


The electrons of energy 105 MeV. have been used for the experimental. 
investigation. The scattering measurements were carried out as described 
in (II.) on a prime cell-size of 254. The readings which were studied 
were those‘on 100, cells with maximum degree of overlap corresponding 
to A=4. 

The tracks were divided into segments of 1400. Eighty-five segments 
were thus obtained, each containing 49 values of D,(A) for A=4. Each. 
segment was regarded as being an independent sample from a universe 
of electron tracks of energy 105 MeV. From each segment two further 
sub-samples were obtained. One consisted of the 13 values of | 
D, corresponding to A=1. The other was made up of 7 values of 
D, obtained by choosing from the samples of 49, (A=4), every Sth reading 


1054 C. O'Ceallaigh et al. on the Observations on the Multiple 


starting with the first. Since the correlation between the individual 
values of D,; in the samples of 7 was negligible, it was expected 
that the standard deviation of their sampling distribution would agree 
-with that predicted by the statistical theory of small samples. It is shown 
later that this was indeed the case. 

To examine the effect of increasing A on the spread of the sampling 
distribution, the statistic chosen was m, ((3), above). Because of the 
small number of signals in the samples, inclusion of the occasional large 
value of D, would have a serious influence on the measured spread. 
‘To reduce this, all tracks which contained as components of the samples 
of 49, one or more values of D, greater than 4 D, were temporarily excluded 
‘from the statistics. It was felt that this procedure rendered more exact 
the comparison between the standard deviations of the different samples. 
‘The results are summarized in Table Il. The symbol E(m,) represents 
the expectation of m,, and is the mean value of the sample values of 2. 
It is the best estimate of , the mean deviation of the universe. 
J», is the experimental standard deviation. The uncertainties quoted 
are the estimated S.D. of o,,.. The upper histogram in fig. 1 shows the 
experimental distribution for the samples of 7 (virtual independence), 
and the lower, that for the samples of 49, (A=4). 


TaB_e. IIT. 
Fifty-seven Samples of Electron Track L=1400,. 


m\ N Kon) Ting 
4 49 0-516 0-099 +0-009 
I 13 0-531 0-114+0-010 


Ee 7 0-525 0-130 40-012 


These results appear quite consistent with those predicted by the 
‘theory of Moliére, when the statistical weight of the experiment is taken 
into consideration. Though the theoretical argument has been based 
on a different statistic, viz., s, it seems reasonable to suppose that the 
validity of the conclusions will be little affected, even for samples as 
small as 7. The results given in Table V. lend support to this 
supposition. The difference between the results for o,, was tested for 
significance by the F-test. (Hoel, Introduction to Mathematical 
Statistics, p. 152.) Considering the results for the independent samples 
of 7, it was found that they differed significantly (at a level of ~3 per cent) 
from those for the samples of 49, (A=4), but not from those obtained for 
the samples of 13, (A=1). Bearing in mind the trend of the results and the 
appearance of the histograms, we conclude that a small but definite 
sharpening of the sampling distribution occurs as A is increased, but 
that the degree of improvement which results, is quite out of proportion 
to the increased labour. These experimental conclusions agree with those 
to be drawn from the theoretical treatment of Moliére. 
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Upper histogram. 


2:0 30 4:0 a @) 6-0 70 8:0 90 
Lower histogram. 
Experimental distribution of the statistic mj (mean deviation from assumed 


universal mean=0) for 105 MeV. positrons, 100u cells. Lower histogram 
for samples of 49, (A=4). Upper histogram for samples of 7 (independent). 
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§5. SAMPLING DISTRIBUTION OF SAMPLES OF 7.—CoOMPARISON 
witH THEORY OF SMALL SAMPLES. 


Our samples of 7 were so chosen that they were composed of 
independent values of D;. The sampling distributions of the statistics 
should agree with those predicted by the theory of small samples. An 
experiment on the decay spectrum of »-mesons at present in progress, 
depends for its interpretation on a knowledge of the form of the sampling 
distributions. (O’Ceallaigh 1951). It is of interest, therefore, to test 
experimentally the correctness of the above assumptions. 

Eighty samples of 7 were chosen. They differed in one respect 
-only, from the 57 samples already discussed. As before the signals 
D,;>4D were rejected, but the loss was now made up by including the 
value of D,; next following. The values of all three statistics defined 
in (3), (4) and (5) were computed for each sample and their distributions 
were plotted as histograms. The standard deviation of each was 
calculated. 

The observed distributions are reproduced in fig. 2. The upper 
histogram which gives the distribution of the standard deviation s, (5), is 
compared with the theoretical distribution of the 8.D. of samples 
of 7 drawn from a normal universe. This is given by the well known 
expression of (Kenney 1940, Mathematics of Statistics —IL.) : 


2)3(n—1) 
his)= Oe (—4ns?/o?)s"—, 1 ee a 
where o is the standard deviation of the universe, and s that of the sample 
of size n. 

_Experimental values of mean and standard deviation for 100, cells 
are derived from the distribution of (I.) fig. 3. These values, respectively 
537 and 0-690 scale divisions, are used throughout the present paper: 
as the best estimates of » and o the universal parameters. 

The lower histograms represent the experimental distribution of the 
mean deviations of m, (4) full line, and mp, (3) broken line. The theoretical 
distribution for m for samples of 7 from a normal universe has been 
calculated from the tables of Godwin (1945). : 

Agreement between the experimental results and the theoretical 
curves is satisfactory. There is some evidence that the histograms 
represent distributions more peaked than would be obtained from a 
normal universe. This does not seem unreasonable when we remember 
that the universe from which the samples were drawn is, in fact, somewhat 
more peaked than the normal. While the position of the mode of the 
distribution of m (full line) agrees well with theory, that for the distribution 
of m, seems to have a higher value which appears to coincide with the 
mathematical expectation of the statistic (cf. Table IV.). 

The values of the standard deviations of the experimental distributions 
are given in Table IV. where they are compared with the theoretical. 
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values. That of the distribution of S was obtained b 


y substituting 
N=7 in (9). 


That for m was found from the expression (Kendall 1947) 
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ex fie distributi isti s of 7 (independent)— 
distribution of the various statistics for 80 samples o 
iota. for 105 MeV. positrons, 100y cells. 
i i istributi tandard deviation with 
arison of the experimental distribution of the s 
ie Seat distribution obtained for samples of 7 drawn from a normal 
universe. ; ah 
i iving the experimental distribution of: the mean deviation, 
q bea fou (fall line), and mo, from assumed universal mean—0 
(broken line). Curve gives the theoretical distribution as calculated from 
the tables of Godwin (1945), for m, for samples of 7 drawn from a normal 
universe. 


SER. 7, VOL. 42, NO. 332.—SEPT. I9Q51 4c 


1058  (. O’Ceallaigh et al. on the Observations on the M ultiple 


The theoretical expected values of the statistics, namely those the mean 
of the individual values of m, m, and s are also given, and compared with 
experiment. 
TABLE IV. 
o=0:690 ; p=0-537. 


Statistic S8.D. Sampling Distribution | Mathematical Expectation 


Experiment Theory Experiment Theory 
Mog 0-153 0-541 0-537 
m 0-154 0-156 0-504 0-509 
s 0-178 0-177 0-612 0-613 


Agreement between experiment and theory is better than might be 
expected. We may conclude that when the individual values of D; are 
chosen, as in these samples of 7, they are effectively independent. We 
may safely use the customary theory of small samples to estimate the 
spread of the sampling distributions of statistics derived from them. 


§6. EXPERIMENTAL DISTRIBUTION OF SAMPLES FOR A=1. 


To obtain results of greater statistical weight than those given for 
A=1 in Table III., samples of 7 were selected instead of the samples of 
13 already discussed. They were chosen exactly as were the samples of 
of 7 independent readings, except that every fourth instead of every 
eighth signal in each track was used. There resulted 157 samples of 
7 for which A=1. The distributions of m, and s were plotted, and their 
standard deviations were calculated. The results are given in Table V. 


TABLE V. : 
8.D. of Sampling Distributions of Mean and 8.D. for Samples of 7. 


Condition of Sampling ~—m,y 8 Number of Samples 
Independent 0-154 0-178 80 
A=1 0-162 0-189 157 


If we apply to the results for s the argument given on p. 1053, we find that | 
89 per cent of the total number of second differences obtained with 
A=1, may be regarded as being independent when computing the S.D. 
of the sampling distribution by means of (9) or (10). This figure of 
89 per cent is in exact agreement with that deduced from the theory 
of Moliere. The percentage change in the §.D. of the sampling 
distribution of m, appears to be about the same as that of s. 

We note, however, that the difference between the experimentally 
determined standard deviations is not significant at the 5 per cent level, 
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Comparison between the histograms of the experimental distributions 
for independent cells with those for A=1, seems to indicate that the 
experimental distribution for A=1 is less peaked than that for the 
independent samples, and thus are more in agreement with the theoretical 
distributions for samples of 7 for a normal universe. 


REFERENCES. 


FisHer, 1920, Mon. Not. Roy. Ast. Soc., 80, 758. 

Fow ter, 1950, Phil. Mag., 41, 169. 

GopwIn, 1945, Biometrika, 33, 254. 

Host, 1946, Introduction to Mathematical Statistics, p. 152. 

KENDALL, 1947, The Advanced Theory of Statistics (Griffin), 3rd. Ed. 

Kenney, 1940, Mathematics of Statistics —II. (London: Chapman and Hall), 
(a) p. 136, (b) p. 135. 

MourerReE, 1951, to be published. 

O’CEALLAIGH, 1951, in preparation. 


[ 1060 ] 


CVII. CORRESPONDENCE. 


The 7-Meson. 
By P. EK. Hopeson, 


Imperial College of Science and Technology, London*. 


[Received July 27, 1951.] 


In 1949, Brown, Camerini, Fowler, Muirhead, Powell and Ritson reported 
the discovery in a photographic emulsion exposed to the cosmic radiation 
of a particle having a mass about a thousand times that of the electron. 
On coming to rest, it decayed into three singly charged particles, one of 
which was certainly a 7-meson, and the others probably 7-mesons also. 
A year later, Harding reported two more of these events, which are 
generally known as 7-mesons. They were found in plates exposed in 
a pipe sunk in a glacier on the Jungfraujoch. Recently, Powell (1951)+ 
has given details of a fourth 7-meson. 

The fact that two 7-mesons were observed by Harding in plates 
surrounded by ice, while only a comparable number was observed in 
the vastly greater volume of emulsion examined by the Bristol group 
and others suggested that perhaps r-mesons are more likely to be emitted 
from hydrogen. Those produced in nucleon-nucleon collisions inside 
larger nuclei may have a high probability of absorption before escape. 

A further batch of 400u Ilford G5 plates was therefore exposed in the — 
glacier and one further 7-meson has been observed. A facsimile projection 
drawing of the event is shown in the figure. 

The track of the r-meson dips too steeply in the emulsion for any 
measurements to be made onit. The tracks of the three secondary 
particles, none of which ends in the emulsion, are coplaner within 2°. They 
have ranges in the emulsion of 445u, 720. and 990u, respectively. The 
multiple scattering of the longest track was 1-25-+0-5°/100u. The 
values expected for protons, 7-mesons and p-mesons of the same rate 
of energy loss are 0-2°/100u, 1:4°/100u and 1-9°/100u respectively. It 
is therefore almost certain that the particle which made this track was 
a m-meson or a z-meson. The energies of the particles producing the other 
two tracks were too great for reliable scattering measurement to be made 
on the length of track available. Assuming all the secondary particles 
to be z-mesons, their energies were found to be 33-5-+5, 3144 and 
9+1:5 MeV. from grain counting. Assuming the 7-meson mass to be 
274+2m,, this gives the mass of the 7-meson as 969+15m,. This 
result is in good agreement with those of the Bristol group and Harding. 

If momentum was conserved in the decay process, it was calculated 
from the angles between the directions of emission of the secondary 
particles that their energies were in the ratio 1 : 0-78 : 0-29. The 
observed values are in the ratio 1 : 0:93-0:18 : 0:27+0-06. This 
a gs ee 


*Communicated by Sir George Thomson. 
+ A report of this work appears in this issue (page 1040).—Editor. 
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agreement, together with the fact that the tracks are coplanar, shows 
that there are no more than three particles emitted in the decay process. 
The total energy released is 73-5--7 MeV. 

It is a curious fact that of the 7-mesons reported by Harding, one 
came upwards and the other sideways. The new 7-meson was also 
coming upwards when it entered the emulsion at an angle of 47° with 
the vertical. This can be understood if the r-mesons emitted downwards 
are more energetic than those emitted upwards. They will then take 
longer to come to rest and are more likely to decay in flight than those 
emitted upwards, which come to rest quickly. 


Facsimile projection drawing of an event interpreted as the decay of a 7-meson. 
(Observer : P. E. Hodgson). 

The event described in this letter provides further confirmation of the 
existence of 7-mesons, and also gives some support to the view that these 
particles have a larger probability of being produced in hydrogen than 
in the other elements. Further examination of the plates is in progress. 


- T should like to express my thanks to Professor Sir George Thomson 
for a valuable discussion of this event, to Mr. R. M. Tennent for making 
the multiple scattering measurement on the track of one of the secondary 
particles, and to Mr. A. J. Herz for developing the plates. 
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Paramagnetic Resonance in Gadolinium Ethylsulphate. 


By B. Buraney, F.R.S., R. J. ELLiorr, H. E. D. Scovin and 
R. S. TRENAM, 
The Clarendon Laboratory, Oxford*. 


[Received July llth, 1951.] 


Kereiaar’s (1937) determination of the crystal structure of the rare 
earth ethyl sulphates indicates that the rare earth ion is in an electric 
field of Dy, symmetry. On this basis Elliott and Stevens (1951) have 
succeeded in explaining a number of paramagnetic resonance results. 
In the case of Gd+++ the ion is in an 8S state, which is affected by the 
crystal field only through high-order interactions which are difficult to 
evaluate. We may assume, however, that the energy levels should 
fit a spin Hamiltonian of the form 

H =gBH.S + AP2(8)+ AYPy(S)+ AgPe(S)+AgPe(8)  . (1) 
which reflects hexagonal symmetry about the z-axis. Here P7"(S) is a 


form of operator which has the same transformation properties as the 
corresponding real surface harmonic Y;"; typical forms are 
P§=38?7—S(S+ 1), 
Pp=i(84. +82). 
The coefficients Aj” are related to the crystal field but not in any simple way. 

Previous work on gadolinium salts has always been interpreted in 
terms of a crystal field of predominantly cubic symmetry. We have 
therefore carried out a number of measurements on gadolinium ethyl- 
sulphate diluted with lanthanum ethyl sulphate to verify the correctness 
of the Hamiltonian (1). Experiments at 1-3 cm. wavelength showed that 
the spectrum consists of seven roughly equally spaced lines, whose 
. splitting varies approximately as (3 cos? @—1) with the angle @ between 
the external magnetic field and the crystalline hexagonal axis. This 
shows that in the Hamiltonian the predominant term is PS’, with axial 
symmetry, not cubic symmetry. Thus gadolinium is comparable with 
manganese, where the splitting of the °S state, previously interpreted by a 
cubic field, is known to be mainly axial in symmetry (Bleaney and 
Ingram 1951). 

To determine the values of the higher terms in (1), measurements were 
inade at 3 cm. wavelength in an apparatus incorporating a proton reson- 
ance system for determining the magnetic field accurately. From 
measurements at 20° K. with 6=0 (magnetic field parallel to hexagonal 
axis), the constants are determined as (in units of 10-4 em.~!) 

3A2=+200:3 ; 60A9=—3-93; 1260A°=—+0-52. 
The lines are symmetrically placed about the central line within 


2x 10°4 cm." showing that off-diagonal elements of the Hamiltonian 
(such as P§) are small. 


* Communicated by B. Bleaney, F.R.S. 
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The spectrum observed perpendicular to the axis (@=90°) is markedly 
asymmetrical. This is due to second order effects which agree within a 
few units of 10-4 cm.-! with those calculated from (1) transformed to 
take the new direction of the field as the axis of quantization. These 
second order effects can be eliminated by taking the mean of the corre- 
sponding intervals on the low and high field sides of the spectrum, and a set 
of constants can be determined independently of the previous measure- 
ment except for a small third order correction. These are (in 10-4 em.~!) 


3AS=+201-3; 60A%=—3-77; 1260(2A9—1A%)=+0-41, 


- Comparison with the previous set shows good agreement for A? and A°, 
and we are justified in using the previous result for A? to obtain 1260 
Aj=—40. The value of gy has been determined in two ways: (a) 
comparison of the $<—>—} transition with the resonance line of diphenyl 
tri-nitro-phenyl-hydrazil (cf. Bleaney and Ingram (1951)) at 1-3.em. 
wavelength, which gives y=1-992 +.0-002 both parallel and perpendicular 
to the axis; (b) comparison with the proton resonance at 3 cm. wave- 
length which gives g (parallel)=1-:994+0-001, g (perpendicular) 
=1-994+0-002. The error is greater in the latter case owing to the 
second order effects which are more important at 3 cm. and which may 
be due in part to small off-diagonal terms omitted in (1). 

In an effort to elucidate such terms the lines have been plotted into 
zero magnetic field at long wavelengths. Owing to the Kramers 
degeneracy they reduce to three, at 0-1127, 0-0834 and 0-0489 cm.-! ; 
the values calculated from the proton measurements using (1) are 0-1130, 

-0-0833 and 0:0455 cm.~! respectively. The discrepancy in the last is 
well outside the experimental error, but so far it seems impossible to 
account for it by adding terms to (1) (such as P?, Pj, which reflect trigonal 
symmetry) without upsetting the agreement with the spectrum observed 
at the shorter wavelengths. 

The reduced population of the upper levels at 20° K. in the strong 
fields necessary to observe resonance at 1:3-cm. wavelength enables us 
to determine the signs of the coefficients as given above from the asym- 
metry in the intensities of the lines. Hence at zero field the levels lie 
at 0, 0-049, 0-132 and 0-245 cm.-1, the +7/2 levels being uppermost. 

Measurements at higher temperatures show that the overall splitting 
is about 11 per cent less at room temperature than at 20° K. At 90° K., 
the gadolinium resonance has been observed as an impurity in crystals 
of the ethyl sulphates of cerium, neodymium and samarium. ihe 
splitting decreases as the atomic number of the predominant ion increases, 
the change being equivalent roughly to a decrease of 3 per cent in the 
coefficients A, A? per unit increase in atomic number. 

Preliminary measurements on gadolinium bromate, nitrate and sulphate 
show that the spectrum is of the same nature as in the ethyl sulphate. 
The splitting is rather smaller in the bromate, and larger in the nitrate 
and sulphate than in the ethylsulphate. 
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As previously reported (Bleaney and Scovil 1950) no hyperfine structure 
due to the odd isotopes of Gd could be observed, and it must therefore 
be obscured by the tails of the strong lines due to the even isotopes. 
In view of the large hyperfine structure found in Nd and Er, this indicates 
that the 8S state of Gd must be substantially a pure 4f7 state. 
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In this period when physics is advancing so rapidly, and when so many contri- 
butions to original research appear in the scientific journals, there is a very great 
need for summarizing articles which will allow scientists, without too great an 
expenditure of labour, to find out what is going on in subjects akin to their own. 
Physicists are therefore indebted to the Physical Society for continuing to 
maintain the high standards of Reports on Progress in Physics, which is at 
present the only journal in this country which contains such summarizing reports. 
The present volume contains eleven contributions ; these deal with spectro- 
scopy, including optics, with the physics of the atmosphere, and with nuclear 
physics and cosmic rays. The volume ends with an article by Professor 
Mayneord on applications of nuclear physics to medicine. N. F. M. 
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